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FOREWORD 


This report was prepared under contract NAS 8-114-95 and is one of a series 
intended to illustrate analytical methods used in the fields of Guidance, 

Flight Mechanics, and Trajectory Optimization. Derivations, mechanizations 
and recommended procedures are given. Below is a complete list of the reports 
in the series . 
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1.0 STATEMENT OF THE PROBLEM 


This monograph will develop near optimum steering equations for boost 
vehicles. However, a very severe constraint will be imposed which will 
remove this problem from the spectrum of problems normally formulated in 
a rigorous manner and solved numerically on a large digital computer (see 
monographs on optimization formulation techniques, SID 65-1200-4 and SID 
65-1200-7). This constraint is that the solution must be possible in less 
than real time* on a small digital computer of the type employed in the 
guidance system of boost vehicles. 

Three separate approaches to boost guidance will be developed and 
discussed. These are: 

1) Path Adaptive (iterative) Guidance 

2 ) Explicit Guidance Employing Guidance Polynomials 

3) Perturbation Guidance 

These discussions will have as objectives the demonstration of differences 
in the formulation of the guidance problem, the development of a means of 
assessing the mechanization requirements** for each, and the presentation of 
information necessary to assess the potential accuracy, the flexibility 
and the limitations of these three forms of guidance. 

Path Adaptive (iterative) Guidance as it will be discussed refers to 
the use of an approximate solution to the equations of motion to iteratively 
(that is, the solution will be repeated at points along the resultant tra- 
jectory) define the optimum steering logic. The assumptions necessary to 
generate the approximate solution in regards to the nature of the thrust as 
a function of time, algebraic and trigonometric approximations, gravitational 
approximations, and time-to-go approximations will all be discussed. The 
attention will then shift to the mechanization of the resultant material to 
demonstrate the performance of the technique and indicate sections of the 
logic which could be improved to produce more optimum performance. Path 
Adaptive Guidance has been discussed in the open literature by several authors 
(References 1.1, 1.2, 1 . 3 , 1.4, 1.5, 1.6 and others). However, an approach 
will be taken here which parallels that of two, I. E. Smith and G. W. Cherry 
(References 1.1, 1.2, 1.3, 1.4). Minor differences will, of course, be 
apparent if the material is compared to the references; however, reasons 
for the differences will normally be given in way of justification. 

The second section of the monograph presents a discussion of another 
approach to the development of a path adaptive guidance system. The objec- 
tive of this section is the formulation of the guidance problem in terms 


*Real time is time as observed by a stationary clock. The requirement that 
the solution be possible in less than real time assures that subsequent 
time will be available for corrective action if it is required. 

**No attempt will be made to assure that the same guidance computer can be 
employed in each of the three approaches. Rather, the requirements will be 
used to indicate the type of guidance computer which is acceptable. 
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of polynomials involving the instantaneous state of the system. The major 
motivation behind this approach is the fact that a steering logic which is 
somewhat more efficient than that mentioned in the previous paragraph could 
be obtained at the expense of some additional preflight computations. However, 
other motivation and objectives and techniques employed in formulating the 
guidance equation will be presented in the development. 

The perturbation guidance discussions differ from those mentioned in the 
previous paragraphs in that this type of guidance tacitly assumes that the 
mission is completely defined before launch, that a reference nominal tra- 
jectory is available, and that deviations from the reference trajectory will 
be small so that linear theory can be employed to generate the steering 
commands to return the vehicle to the nominal trajectory. This approach is 
extremely simple since most of the difficult computations can be made before 
launch on a large ground based computer. Further, the resultant trajectory 
will be more nearly optimum than those yielded in the other approaches to 
guidance since the zeroth order steering data (provided on tape) will be the 
result of a simulation of the optimum boost problem on an "exact" model of the 
earth and vehicle. The method has one major disadvantage, however, in that it 
requires large amounts of precomputed data. These data, in effect, limit the 
flexibility of the system once it is launched. 

The discussions of perturbation guidance are divided into two principal 
parts : 

1) Velocity-to-be-gained (or required velocity) 

2) Linearized Perturbation Guidance 

In the first of these approaches, the velocity required to assure a position 
constraint on the terminal state is the controlled parameter. These discussions 
present the development of quantities to be utilized as error signals and of 
the optimal use of these data to generate the required steering commands. The 
works of Sarture (Reference 3.3) and others (References 3-4 and 3*5) are re- 
viewed to define the complete analytic framework of the problem. In the second 
approach, methods for employing instantaneous measurements of the state de- 
viation to generate guidance commands intended to minimize some measure of the 
terminal error are discussed. This portion of the monograph will concern 
itself with the material of the type presented in References 3.19 and 3.20. 

The first and third of these approaches to the guidance problem have been, 
or are being, successfully applied to boost vehicles. The second, due princi- 
pally to numerical problems implicit in the preflight simulation and con- 
struction of the coefficients of the polynomial is not presently being 
considered for application. Thus, major emphasis will be placed on the 
iterative and perturbation guidance approaches. These discussions will present 
detailed derivations of the equations which will be mechanized and will attempt 
to establish the basis for all of the assumptions which have been made. 
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2.0 STATE-OF-THE-ART 


2.1 PATH ADAPTIVE GUIDANCE (ITERATIVE) 

2.1.1 Introduction 


In the other sections of this monograph, a series of methods are dis- 
cussed which are capable of generating steering commands as a function of 
time and the state at that time. These methods were developed upon the 
assumption that the computational capabilities of the guidance computer are 
very limited and thus employ large quantities of precomputed data on the one 
hand or guidance polynomials which are difficult to compute on the other. 

This section presents an approach to the problem which is capable of over- 
coming the major objections to such guidance schemes and which is still 
sufficiently simple as to avoid the tremendous computational requirements of 
a rigorous variational calculus formulation. 

The mechanism of this formulation is an analytic solution to the optimum 
programming problem for the special case of a constant flow rate under the 
assumption that an average gravity vector can be defined. While the results 
of this process are reasonable, provision has been made to allow for the re- 
petitive solution of this problem using the progression of measured states to 
achieve a high degree of accuracy, in the large without requiring each of the 
independent computations to be precise. 

The method of analysis is similar to that first prepared by I. E. Smith 
in Reference 1J; however, minor modifications have been made in the procedures 
to effect a higher degree of similarity between the results of this analysis 
and the results of a precise optimum transfer problem. To aid in the appre- 
ciation of the assumptions employed and of limitations imposed on the analy- 
sis, frequent reference will be made to relate approaches to this guidance 
concept. The purpose of these references will be to contrast the accuracy 
and mechanization differences as well as the relative efficiencies of the 
approaches being considered. 

2.1.2 Out-of-Plane Guidance 


The general boost-to-orbit guidance problem is generally preceded by one 
in which the launch time is selected such that the displacement from the 
desired plane (i.e., the plane containing the terminal state) is reasonably 
small. This step is accomplished in order to minimize the amount of work 
expended to turn the plane of the trajectory at the time the transfer and 
target planes intersect. Thus, the differential equation for the motion of a 
vehicle in this perpendicular (to the target plane) direction can be separated 
from the remaining equations with small error provided a first-order correction 
is applied. This observation is illustrated in the sketches and formulation 
which folio;/. 
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First, consider the definition of the coordinate system 



a /\ A . A 

where R, S, W are unit vectors with origin at the aim point (R along the 
desired radius vector; W along the desired angular momentum vector); second, 
define two orientation angles for the thrust vector in this coordinate system 



where F is the applied force vector per unit mass whose magnitude is F and 
whose orientation is defined with respect to the target plane by the angle 
O' . (This choice of reference axis was made to separate the out-of-plane 
dynamics to as great a degree as possible from the in-plane dynamics. This 
selection was based upon the consideration of the fact that (X (t) will be 
quite small for most trajectories since the plane of motion can be controlled 
to a degree by launch time selection and since the terminal value will in most 
cases be very nearly zero in order to assure that the terminal velocity per- 
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pendiculaf to the target plane will be zero. Under these definitions, the 
equations of motion normal to the target plane are 

W = r-w (2.1.1) 

Now, the thrust attitude angle, , will be defined such that the change in 
velocity parallel to the desired plane resulting from out of the plane thrust 
is as small as possible; i.e., the loss function 

J ~ F (/ - Ctrl &C) dt ( 2 . 1 . 2 ) 

is as small as possible. 

In contrast to a solution of these equations which strives for maximum 
accuracy, the guidance problem in general strives for both relatively high 
accuracy and computational simplicity. Thus, the first of a series of 
approximations will be introduced to allow analytic solution of the problem 
by uncoupling the out-of-plane dynamics. It is assumed that the integral of 
the corrective term in the equation of motion (that containing gravity) can be 
expressed as 



r 


/\ 

vV dt 



G. 


( 2 . 1 . 3 ) 


Attempts to evaluate this time average acceleration will be postponed until 
the present discussions can be concluded. It is noted, however, that since 
the integral itself is small, relatively crude approximations (for example, 
truncated series) for both w and r as functions of time could allow the dot 
product g . w to be evaluated independent of the analysis of the in-plane 
motion with adequate precision. This procedure doesn't appear desirable 
though, due to the fact that subsequent steps (in the analysis of in-plane 
motion) will demonstrate the necessity of generating an accurate estimate of 
the time average gravitational acceleration^ This fact means that G can be 
estimated from the equation G = ga V g . W with more than adequate pre- 
cision. 

Now, since an average value of the gravity correction, G, is assumed 
known, the optimum thrust coxomands for oc — (t) can be determined by 

considering the optimization of 

J(o0 = f f(/ - C&lcK) dt 
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subject to the integral constraints 


j) (u) --w(r) + w(o) -f\A/(o)T-+ [fF^ads^dt (2 ' 1 ‘ 4) 


J Z (<*)*-\A/(T) ■+• G7 + W(0)+[f ^u>t cSe/t ( 2 . 1 . 5 ) 

y o 

But this statement is identical to that of the class of isoperimetric problems 
of the calculus of variations; i.e., if the augmented integral I* is defined 
to be 

I* = 1(d) + \,J, ~h \ Z J Z 


"fs* dt 


then, the Euler-Lagrange equations for the system become 


a/' d_r bj_ 7 '\ __ 
dt \ bck ' ~ 


d_r bfj 
dt \ bck 


But, I* (minus the terms which are not functions of a ) is 

, T T 

I * F(/~ C#d <7\ )t A ,( £ C\ t/jJ dt / A z J T f 'dcA a\ dt (2.1.6) 


This equation can be simplified by integrating the second term by parts as 

I J c< dsdt - T f o\ eft - ft cK c/t 

o S 0 S n 


This fact allows I* to be written as 

,T 
'O 


I*=f [f ( / - Ctd <x) + ( A,T + A Z )(FS^(^)~ X,tF S^c/6 ( 2 . 1 . 7 ) 


v or f* as 


f* = F( /'£**<*)+ Fh,(T-t) + A Z J SLdrucK . (2.1.8) 
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Thus, since f* is independent of « , the Euler- Lagrange equations are 







and since F ^ 0, 


25*-*- CK 


- A 2 

- ( A /7 y- A 2 ) + A, t . 


or 


= A y Bt 


,4*orL cA 


A + Bt 
IS/tCA-tst) 2 


(2.1.7) 


independent of the nature of F. 

The values of A and B can be obtained by integrating the equations of 
motion (i.e., the constraint equations) and substituting the boundary con- 
ditions for the terminal and initial states. This task can be performed for 
the general case, however, the resulting expressions are nonlinear in the 
constants A and B. This fact means that relatively elaborate procedures are 
required to evaluate the constants and suggests that a simpler solution would 
generally be preferred for the on-board guidance problem. Thus, an attempt 
will be made to approximate sin oi for the case where the change in oC is 
relatively small (relative to one radian) during the burn of any single stage. 
For this case 

t&sfC o\ - (A* t ACk) 


and 


A ck 


A / Bt ~ CY 
^ a Ok* 


(2.1.10) 


where for the moment o< * is an unknown constant reference. Now since Ack. 
is known, the sine of «. can be written as 

s&tsr'L (fa = C Ok* -A A Ck ) 


C** + A A4K1 Ok 
CF* A 

C ¥ 


( a y Bt - ± 


( 2 . 1 . 11 ) 


2 cX * 
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which again is linear in time. This fact means that a linear sine steering 
program will be near optimum provided that the assumption (A«=s) z << I is 

not violated. For this reason, the degree-of-optimi z ation will be com- 
promized slightly to effect major simplifications in the guidance process and 
the constants C and D will be defined by fitting the boundary conditions to 
the equations of motion. Note that the angle « * which was introduced for 
the purpose of expansion is not required nor can it be from this information 
(2 equations and 3 unknowns). 


Consider the case where the motion of the vehicle is produced by a con- 
stant flow rate motor (this type of operation is characteristic of nearly all 
large liquid motors and most large solids); i.e.. 


F 


F( 6) E Fo 

I ~ nr\(_t +tb) (/-/yif" ^ t 


( 2 . 1 . 12 ) 


where t b denotes the burning time which has ellapsed since this particular 
stage was ignited and where both F and m have been normalized by the initial 
mass of the stage /i.e., F = thrust/m(o); m=mass rate/m(o/7 Under this 
assumption, the linear sine steering program will produce the following re- 
sults for the constraint equations. The terminal velocity is 


vO(T) 


FeCC tOt) 


J 1C 1 ~ mi 

- w(o) 4 GT 4 F 0 C 


+ G | db + w[o) 

c Ft 


r<- 






tdt 


(/- 


— (2.1.13) 


= W(0) t 6T t F 0 c{ + F o 0f 2 


where 




r 

nn L 


T 

/- rht b 


ft 


fb /r\ 2 . l~ 


- T 

rb 


+ 


rb 
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And, the terminal position is 
(j) - f T f t \ F ° ° - S ~- + & \dSdt + W(o)T + IA/CO) 

L ' J L Vt-AW-AS J 


O o 


- \JJ{ 0 ) t wCo) 7 i- - T 2 + F 0 c 
2 


T t 

III 


c/s dt 


* 


T.t / _,y (2.1.14) 

j c/s dt 


<i-*ib)-ri<s „ 0 


J f “** 




- y ui/(0) / / — T 2 t F a C G t 

where 2 * 


+ F ° 


J, -f^OM -/- + **{' -T^] di 


/ '"Mf 

mT "I 

jtn ^ 

/- ~ r ) 

- / 

t / 

m l rn /I 

J — . 

!-rht b \ 

/ - mt b y 




3* 


f-(t)clt 


.m 

y L 


, /, I dt 


- 

c 

I 

. 


T‘ 


/ ~ nn t, 


2m 


rr\ 


Thus, since these equations are linear in the unknowns C- and D, the solution 
to the system 


wCt) ~ W(o) - 6 7 *] 

1 G7 Z I 

►- F a 

f, 

fz' 

- 

c ' 

k WL7) ~ W(o) - \VC°)7 " J 


V 

2* 


o 

i 


can be evaluated as follows: 


c 

1 

' . -tffc 

-^1 

WCT) 

* 

- w(o) 

~ G 7 _ 

. „ , GrT* 

D 

v J 

** b, 

#■- J 

WCT) 

- w Co) 

- W( 0 )T - — 

2 > 


where 


( 2 . 1 . 15 ) 


4 = f,<3 2 -f z %, 
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There is one chance for error in this solution resulting from the fact 
that a -*■ 0 as T -*■ 0. This fact requires that the limits for both C 
and D be established in order to assure that the solution is always deter- 
minate. This step is not readily accomplished. However, by noting that in 
the limit W, W, W and G approach zero (the terminal plane is the desired plane), 
it is possible to conclude that C and D should also approach zero in the limit. 


This solution completely solves the problem under the assumption that 
sinoC — C+Dt. Further, though there is no assurance that the approximations 
which allowed this solution will also be reasonable, it is noted that suc- 
cessive re-evaluations of C and D as the epoch corresponding to zero time is 
moved along the trajectory (in terms of w and w) will assure that the process 
will converge for all of the cases of interest. However, the degree-of-opti- 
mality for cases involving values for which the sine of oe. cannot be approxi- 
mated as specified is questionable) . 


A slight variation of this approach was employed by Cherry (Reference 
1.3) 1.4). However, rather than assuming the nature of r(t) for the 
purpose of computing the gravity correction. Cherry noted that the linear 
tangent steering which results for G = 0 or a constant is approximately 
equivalent to linear sine steering. He then assumed a "near optimum" steering 
law 


& X 

AtX <h - C+Pt / 


Atr tv 
r 3 


¥ 


This assumption causes the gravity contribution to be cancelled from the 
equation of motion and allows C* and Eh' 5 - to be evaluated in a manner similar 
to that outlined on the previous pages. However, if a reasonable estimate 
of G can be made, the approach outlined previously will always provide a 
better estimate of the optimal <X = (t) with only slight additional 

complexity. 

2.1.3 In-Plane Guidance 


The discussions of the previous section provided the basis for analyzing 
the motion normal to the desired target plane (the plane of r^, ’v-f). Thus, at 
this point the analysis will be restricted to the motion in, or parallel to, 
the desired plane. Consider the sketch and nomenclature presented below: 
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SJL A /\ 7 

The equation of motion in the instantaneous coordinate frame [r, W x. vj 

5 u 


S = - + F Coo C* Coo 

r 


(2.1.16a) 


U - 



Coo C*- /S 


(2.1.16b) 


where 



r 


(2.1.16c) 
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r= U 


(2.1.l6d) 


and where the approximate form of the optimum program for «■ (m = constant) 
was derived in the previous section. 

If these equations /denoted by f/x(t), u(t), t / - x = Oj are to be inte- 

grated along a path which yields a minimum value for some scalar measure of 
performance (for example, propellant expended) defined by the equation 


j = <p\x(T ) , r] +J T r[z(i) i u(t),t]dt ( 2 . 1 . 17 ) 

4 

where u (t) represents the control vector, then the differential equations 
defining the optimal path (control) can be generated by adjoining the system 
of differential equations to the performance index 


J* <P\*(T),T] *A r r xU ) , u(t) ,t] + i r (t){f[zU), u U),t] 

- z J ] dt 

= <p\_ %(T),r] + j U-/-\x)d£ ( 2 . 1 . 18 ) 

4 

where H denotes the Hamiltonian. But the integral of this equation can be 
evaluated by parts 


I 


(H-'Ax)dt 




so that 


J = <p [*cr) t l]+ \ t (T) x (T) -A r (tJxU 0 ) 

[x(t) , u. Ct) , t] + £(t)z(i)Jdt 


( 2 . 1 . 19 ) 


Now J is to attain an extremum along the path which is produced by a special 
control u = u(t). This condition is obtained by forming the first variation 
of J. 


= [(!x ■ A sx ]/ kT£ *I + / [(^ + A ' T ^ X4 i^^J di 

v* ft t- 
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r 


At this point in the solution, the multipliers X are arbitrary. Thus, they 
can be selected in such a manner as to simplify the problem. Note in the 
previous equation that variations in the control and in the present state 
can be specified arbitrarily. However, the variations in the state (.X) at 
times in the interval /t D , T 7 and at T are dependent /i. e. , no a’ priori 
information is known. <5 X(tj/ . Thus, the multipliers will be selected such 
that the coefficient of & X(t) is zero for all t 0 - t - T. 

r — 


with boundary conditions 


dX 

1L- r df 

dx d7C 

(to simplify the bracketed term outside the 


( 2 . 1 . 20 ) 


integral) 


i[ . (T) - 


d(p 

d *i 


( 2 . 1 . 21 ) 


for all unconstrained components of the state at t = T and 6X (T) = 0 for 
the constrained components of the state at t = T. Under these conditions, 

<5 J reduces to , 


Thus, if 


<5J =• A Ct 0 )6xCt 0 ) 



dH 

da 




<5 J = 0 for arbitrary 6X (t Q ) and 6u(t), then 


( 2 . 1 . 22 ) 


dH 

9U 


O 


(2.1.23) 


for all t 0 t t - T. 


Now consider the special case of interest in the guidance problem where 
minimization of fuel expended for the case of a constant flow rate is the 
objective. For this case, the problem reduces to the brachistochrone problem 
(i.e., the comparison function is simply the terminal time). Thus, the 
Hamiltonian is a function of the end conditions alone. This fact allows H 
to be written as 


H - A ( ( F coo ot coo /3 - 


SUL 

r 


) 



Fcoo (X Avn. fi 



+■ 



+ \(U) 
i 


(2.1.24) 


Now, since F and cx are known functions of time (for assumptions made 
earlier), the only control possible exists due to the angle /5 . Differ- 

entiating the Hamiltonian to obtain the equations defining the optimal policy 
for /3 yields 


— -O = - X, F coo a. Mm. jS + 'X, FcoocX coo A 
d/S ' * 
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or 




where X 2 and X ^ are defined by integrating 


U ^ ZS I 

X, — X — ■ X 

i i r A r 3 r 


( 2 . 1 . 25 ) 


(2.1.26a) 



(2.1.26b) 



(2.1.26c) 


. SL L 

\ = r z _X x 




(2.1.26d) 


The problem encountered at this point is one of computational complexity 
in spite of assumptions made separating the out-of-plane and planar dynamics. 
This is not to say that solution is impossible; rather, solution seems un- 
necessarily complex for application to a small on-board guidance computer. 

Thus, attention will turn to further simplification of these equations which 
will allow an analytic solution to be developed. Subsequent discussions will 
then attempt to assure that the accuracy obtained from a mechanization of the 
analytic solution will approach that obtained from a more rigorous formulation. 

The first observation is that the range angles subtended during powered 
flight for systems in which the thrust acceleration is of the order of the 
gravitational acceleration are generally small. This fact means that the 
direction-of-gravity is nearly constant and suggests that a "flat earth" model 
can be employed to advantage. Consider the following sketch showing the 
initial position vector, the target position vector and the set of unit 
vectors (R, S) which could be utilized for the purpose of constructing a flat 
earth model (this set conforms to the notation of the out-of-plane guidance 
discussion) . 
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r 



In this coordinate system (R, S), the direction-of-gravity is unknown as a 
-function of time. However, as was the case in the analysis of the out-of-plane 
motion, an average gravitational acceleration can be assumed. Unfortunately, 
the solution of the problem is relatively sensitive to errors in this assumption. 
Thus, care must be exercised to assure that both the magnitude and direction of 
this acceleration are selected in a rational manner. The estimate recommended 
for the generation of this model is, therefore, the time average. Subsequent 
discussions will define the procedure to be followed in the generation of this 
time average gravity vector. Thus, the present discussion can assume that the 
time average has been developed and is known. 


The first step in the derivation of analytic steering equations is now 
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Under these assumptions, the approximate equations of motion are. 


- F CcV CK $ (2.1.27a) 

V r = (2.1.27b) 

s = ]/ s (2.1.27c) 

r = |/ (2.1.27d) 


So that the Hamiltonian for this simple case is 


H = t^-Js)i\ v ^Fu^cA^/ 3 /f ti ^F)+X s (y s ) f A r Lv R ) (2.1.28) 

tf 


and the differential equations for the time varying Lagrange multipliers 
(obtained from 


T - d H 


are 



A\ 


0 



o 


But these equations are uncoupled and may thus be integrated directly to 
yield 


X 


r 




-C,i *c 3 




So that the form of the steering program discussed earlier becomes 


A v _ -Cgt + Cg 
txvn A ~ : — “ „ 

A., -c,t I-C 3 


v s 


(2.1.29) 
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(This law is commonly referred to as the bilinear steering law.) The four 
constants of this law must be evaluated from the boundary conditions imposed 
on the problem. Before performing this evaluation, however, it is recalled 
that the boundary conditions for the can be selected arbitrarily if the 
terminal state is not constrained. [_ (This fact affords the possibility of 
simplifying the form of the steering law without attempting to solve the 
equations of motion.}/ For example, 

a t (T) 


Thus, if 

<P*<P[ scr;J , x s = c,--o 

Under this condition, the bilinear steering law reduces to the linear tangent 
steering law. Note that relaxing only one of the constraints on the terminal 
state eliminated two of the four constants of integration (the two new con- 
stants appear as ratios). This fact and observations of the small sensitivity 
of the resultant traiectories to the constraints on the longitudinal position 
component (noted in References 1.1, 1.2) have lead several investigators to the 
assumption that the performance afforded by the bilinear law can be adequately 
approximated using a linear tangent steering concept. The mathematical 
justification for this assumption will be presented in the following paragraphs. 

The equations of motion 

S-Fcznc\ i ■ S 

r - F u * ck /3 R 

can be made explicit functions of time to aid the integration process pro- 
vided a series of small angle approximations can be made to facilitate 
solution. These approximations will produce a slight inefficiency (in-so-far 
as propellant expenditure is concerned); however, since the launch time is 
selected so that the initial displacement from the desired plane is small, 
since most of the boost trajectories for orbital, lunar and interplanetary- 
flights are- characterized by small values of angle between the velocity and 
local horizontal vectors at injection, and since to attain this state (in a 
near optimum fashion for trajectories composed of an atmospheric arc and a 
guided arc) requires a reasonably small attitude angle when initiating the 
guidance, the assumptions are believed reasonable. 

The first simplification results from the assumptions that the change in 
the attitude of the thrust vector for control of the motion normal to the 
desired plane is small during any single burning process, and that maximum 
amplitude of the angle itself is sufficiently small as to allow the cosine of 
Oc to be represented as 


d <p I 
dX I 







(2.1.30) 


/ - i(C-f Dt) 2 

^ " fr ) 

where the constants C and D are known from the analysis of the out-of-plane 
control problem. 

The second simplification results from the fact that 


/3 - 


~ C-2 1 / C 




-c,t 


/ c 


can be approximated as 




6 -t A/S 
/- A/& 

tasu/S* /- A/S (/ 4 X ) 


(2.1.31) 


- Zttsc /?* / A/S (s<Uc*-/ 3 *) 

provided the change in the thrust attitude during the period of interest (i.e., 
one stage burn) is small. Thus, it is possible to represent a A as 


*0 - - ~ c ^ * I 

P -C,t t C 3 J 

if. 


Cot f c 


(2.1.3 2) 


* 


where 


- C/V 

-cj 


/ c 




- £ 


C, s&X**/} * 
c* - c, t*-*./? * 


q = c 3 ^ uc .*/** 

(Note that the reference angle, /3 *, i s unknown at this point. Subsequent 
steps will, however, define all required information). Now, since a/3 is 
assumed small 

^3 - jo**, 1/3*+ A/S ) 

» /3 * i- A zoo A * 


= ^atn. /5* + 


-C* x t + C* 


(2.1.33) 


- c 


t + c* 


COO /S* 
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and 


coo /3 - coo( <6* + A /3) 

~ aoo /$* - A/5^ao*o x5 * 


- r* 


aoo A* 


Clt -hC 


— /6 * 


( 2 . 1 . 34 ) 


The third and final approximation is now made by noting that since the 
range angle (from initiation of guided flight to injection) is small and since 
most of the trajectories desired are characterized by relatively small values 
of the flight path angle /sin - (? . at injection into orbit, the 

product A/3 sin xS * is much smaller than ax>/3*. Thus, to the 

first order ao°/6 = coo/3 -if-. This observation means that the equations for 

motion in the horizontal direction (S) are nearly independent of the constants 
of the steering program, or conversely that this component of the injection 
position and velocity can be matched to a good degree without enforcing a 
constraint on the program. Finally, since this is the case, the steering 
program for A can be modified to reflect the relaxation of the constraint 
on S(T) without major error. As was demonstrated earlier, this relaxation is 
accomplished by simply equating Cj to zero. The result of this step is now 

* -f- — ^ y -±r I 4*1 £ K, ~fc z t ( 2 . 1 . 35 ) 

L c 3 C J J 


Thus, the problem of defining the constants C 2 /C 3 ", and /2> * reduces 

to any equivalent problem of determining and fL>. 

The composite constants Kl and K.2 can now be evaluated by integrating the 
equation for position in the R direction and substituting the boundary con- 
ditions. This task, while straightforward, is somewhat involved due to the 
number of terms required. For this reason, a shorthand notation will be 
adopted for convenience. 


f, = 


r 

o- tb)-™*- 


f r rhT ' 

A*) / " r~r 


<Z = 


t c/t 


if- rv t b ) ~ d)t rrs 


;2 


\ - r»T 


-(/- rn 6 6 ) |(V- ~ rh 7^j J 


4 . "J**± £ 
m rr\ 
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4-/1 


t*c/t 


_T 2 tt-nnt b )T (/-ntj,) 


'o (/- At b )~At 2 A 


A 2 


4^2 [d, .±1_ 

A 3 L /- ntj. 


ZA 


C /~ ™t b )T t'-*>t b ) Z f, 


m 


/*) 


; 2 


' f -f, Ci)dt - ~~r f 4*i ( / ) d& 


-U— 

/r> l rr 

' - A t b 
A 2 

/ ' At b 


- d\t b 
m 


A 2 


'vKJK- 

"J 


^ T 


9 * ° /V* <*; 


T* 

ZA 

r 


At b 

A 


3 . 


9 s * L 


+ / 


T 3 /-Ax £ 


<£ A 


2 . m 


\z 


ir 2 z i'-?**) 3 

rr j 2 ^ 
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Thus 


r - 


F 0 


(/- F»t b )~ rhi 


[('“ " )-CDt\(K, ~K 2 t) ( 2 . 1 . 


36) 


yields 


&r B r(T) - r(o ) - f a[/j %T^F a j U>/r t f 3 -[( ’ ~^)K Z (2.1.37) 

Q ■Kiri •* 


= rCT)-rCo) - rCo) r-{ $ A „ S ■ $ 7 2 = F 0 {cd« 2 j 3 - [(/ - |*) o\% 

+ ( / - g £')*j 9 ,} 

But these results can be expressed in matrix notation as 


( 2 . 1 . 38 ) 


1 

"n t> 

Q -i ■ 

► - 

, ^ 


i 

t -r" 


/ 

jn f> 
* -» 


> 

2* 

9i_ 



CD 


O 


O 


o- c 4) 


K, 


K, 


CDfi-O- %‘)f s | ~COf 3 
CD ii \ - Co 3* / (' - §^?, 

K- 




(2.1.39) 


' 2 / 


l »2 


fK t Z 


K. 


J V. J 


Thus the solution for the constants is 
? 



A - A/, A Zi A /Z A^i 


(2.1.40) 
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This solution completes the problem of determining the steering equations 
for any single stage of a rocket vehicle under the assumptions and limitations 
outlined in the preceding pages. Before turning attention to the problems of 
generating accurate estimates of the time-to-go and the average gravity (re- 
quired as . a part of the mechanization), however, it is noted that a major 
simplification of the problem can be effected with slight error for most 
cases simply by assuming that since oc (t) is normally small 


604 CK 


/ 


i.e. , 


C ^ D -O 


This substitution will reduce the solution just outlined to one which is 
identical to that for the out-of-plane guidance. This fact can in turn be 
utilized in a guidance mechanization to reduce the storage requirements for 
the computer by allowing the same block of logic to be utilized for two 
applications. 

As a final comment by way of contrast, it is noted that Cherry (References 
1.3, 1.4) approached the problem of in-plane guidance in a completely 
different manner. Once again he noted that a modified steering logic could 
be prepared which would uncouple the equations of motion and allow him to treat 
the problem of in-plane motion as two problems of two degrees-of-freedom. This 
observation permits the same logic used for determining the steering constants 
for the out-of-plane motion to be used for the remaining problems (without 
assuming the flat earth) and assures that the desired terminal state can be 
approximated to the desired level providing that provision is made to control 
the mass flow rate. The guidance law, however, is not as optimal (in regards 
to expended propellant) as that presented earlier because of the fact that no 
mathematical justification (based on optimization) for the approximate steer- 
ing logic can be prepared. Before discussing Cherry's formulation, it is 
noted, however, that his formulation is less complex than that presented 
earlier due to the fact that the gravity estimate has been eliminated (this 
step results in a less optimum solution). The equation employed by Cherry for 
motion along the instantaneous radius is 


where 


r = 


-_M ( V coo Z) . A 

-j- 4* r CJ>0 & /-> 


Z - sOM- 


£ - jiun 


r • v 

rv 

rj 

rF 
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Thus, If 


- M (Vcj° H'f 


is smaller than the thrust acceleration (the term will normally vary between 
— gg and +2g for most orbital, lunar, and interplanetary trajectories of common 
interest depending on the time along the trajectory and the mission), then^a 
perturbed linear sine (Tangent) steering program will be near optimum (pro- 
vided again that land <x z ^i ). The law employed is 


which yields 


^7 (V coo dr) 

AUK jQ - • — - ~ * +■ C r +p r t 

r z r 


r =C r F +D r Ft. 


Thus, the constants C and can be readily evaluated from the present and 
desired radii and radial rates. The fifth component of the terminal state 
is controlled by matching the magnitude of the angular momentum at injection 
by controlling the time-to-go. (This final step will be discussed in sub- 
sequent paragraphs.) 
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2.1.4 Time-to-go 


Implicit in the discussions which precede this section has been the 
assumption that an estimate of the time-to-go was available. Therefore, at 
this point in the process it is essential that several means of generating 
this estimate be considered. 


The first procedure which is of interest was presented in Reference 1.1 
and 1.2 in a slightly modified form. This approach notes that the desired 
terminal velocity can be matched^ as follows: 

‘ * J T 

J o J 

(2.1.41) 






T 

F eft - 


^ ^ 1/0 ~ 3 ax/ 9 

where g a is the time averaged gravity vector (Section 2.1.5) and T is the 
unknown ?§me-to-go. But F is not constant in direction, rather it is being 
turned to effect the steering which is desired. Thus, there is a steering 
loss which will result. This loss must be estimated to assure that a real- 
istic time-to-go is generated. Consider the following vector diagram and 
define the effective thrust as the component of the thrust vector in the 
direction required to match the terminal velocity under the assumption that 
there is no steering ^ 





To simplify the analysis, the assumption is made that the steering loss 
is dominated by the in-plane steering; in general, this assumption is well 
founded; however, the out-of-plane steering losses can be included without 
major revision. Under this assumption, 



f cos c#~ ^*3 


(2.1.42) 


Thus, the first step in the solution process is the definition of the angle 
F *. Consider the equation 
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-AAA 

and the corresponding components in the R, S, W coordinate system 


A'iff cos « = v rf ~ Vf o - 9 ai 


^ - VS 0 - 3 a ^-- 5 ' 


£ 

A 


Thus 


ta.n/2 * = 


Vi - v/ r 


9 « 




-A 

R T 


V -5< 


S<w<j ‘ -S T 


(2.1.43) 


Similarly, the equivalent out-of-plane steering angle ( oc #) can be defined, 
though it is not required at this time, as 


S/A/ cn - 


/ (■ 

L H ' 


v w 0 ~ 


(2.1.44) 


The second step in the process is the expansion of the function 

( fi ~ 0 *) 

o_s (yfl- *3 = ^ f .say * 


where 


= Co 


C05*-/3* - 


S/Z* coS/i + s/JV/3 tasi/3* 


/ 


/ / ta« d # 


(2.1.45) 


-S/A- /3 - K, ~ K z -t 

Cos> /3 - f/ - C K, - Kg t ^ 

~ ' ~ 7 c*, 

= (/-— ) / K,K z t ~ 


L 2 -2 


(2.1.46 


= a. -h t> t -f ct 2 

Finally, the^change in velocity in the direction defined by the velocity 
constraint ( A v eff) can be estimated. 


A V t fl - f F C os C /& ~/& *) c/ 1 £>V e (( 


| A v eff 


' T Fq COS/3* 

(/- /b t 6 ) - sh t 


+ bt + ct 2 + di 
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where 


/ C/~ r>t b )~ rht 


(2.1.47) 


Thus 


where 




(c/+ e.£ + ct Z ) dt 


d = a + Kftctrt / 3 * 

= (' K,tan 0 * 

& - b - d z tan /S^ 

= K , Kz - K 2 tan /& * 
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+ e 


f c 
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/- nhT 

l- dit. 



s Av .0^ - 3 T ~ hT " 


(2.1.4B) 
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Now these two non-linear equations for A V /J2.1.1+1) and (2.1.48)7 
involving the unknown Time-to-go and several functions thereof (g a ver> K^, K^) 
must be solved. This solution can be accomplished iteratively by xxrst noting 
that the steering losses and the gravitational losses, while significant, are 
not dominant. Thus, a first estimate of the time-to-go can be obtained by 

— \ f — - 

- >4)'(v f -v 0 )] 

f mT 1 ^.W) 

IjA / 

L /-«! 


assuming 


i.e. 
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Fo 


rr\ 


J 


T = _ 

~ : / ™ €. 

m L 


_ m ai/. 


«Vf. 


(2.1.50) 


(Subsequent passes through this set of equations for later estimates can use 
the previous solution minus the elapsed time since that solution, i.e.. 


ma.t e 


Tn- t ~ A* ) 


Now, under the assumption that this estimate is sufficiently accurate, the 
iterative nature of the solution can be eliminated by expanding the terms on 
the right hand sides of equations (2.1.41) and (2.1.48) in Taylor series, 
collecting terms to form a polynomial in the correction to the estimated time- 
to-go, and solving. First, define the true time-to-go in terms of the estimate 
and the correction 


T 


= T + 

c. 


AT 


(where A T will in general be small since the corrections are small). Next, 
expand the terms in the respective equations through the second order for 
accuracy 
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and construct the equation 

^ C T „ ) A V.fjr .. C To) ~ *5 (-To) / a ~ h ( / o) 71 z -f £ / o) 


F 0 AT 


( '-'vt-b) ~ ™ to. 


/ 1 •TClo)F_mA/ 2 


gCTo^AT ~ 2 h CT 0 ) T 0 AT + —{ r 

J 2 [ [(!-^b)'^T 0 ] 


- Zk<.T„) 


AT Z • 


V { ~ v o~^T a 


(2.1.51) 


(V f -\/ 0 -^T 0 ) -^1 f -^-] AV^ Z V ° ^ ‘ ] ■ 

J 2 av*(. d AV J 

Now collecting the coefficients of the powers of A T, it is possible to 

write an equation of the form 


A T 

A V L 


A AT 2 t 8 a T y-C - o 


(2.1.52) 


where 
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2 ^ V> 2 A 1/ * 


- h 


The solution is now 
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l^f ' -f ' ) 
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C t- riifi) - rv T 
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'3 ■? /e 2 - vac 

2A 


(2.1.53) 


where the sign ambiguity can be resolved by selecting the root which most 
closely agrees in an absolute sense with the linear estimate of the correction 


^ ' L/ HEAR 
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This procedure will enable the iteration process to be abandoned and will 
generally produce the desired accuracy. Note is made that a process similar 
to the one just outlined was discussed in References 1.1 and 1.2. However, 
the approach presented here is considered superior since uniformity in the 
order of the expansions used in the reference was not maintained, and since 
numerical experimentation has disclosed that many combinations of initial and 
terminal conditions require significant steering (i.e., the steering losses 
are significant so that the time-to-go estimates while reasonable are not 
accurate). Though it is noted that the desired terminal state can be achieved 
in spite of the resulting inaccuracy by the simple expedient of holding the 
steering angles at their maximum value (+) until the estimated time-to-go 
becomes reasonably accurate, the "fix” procedure will result in a much less 
optimum steering program. 

At this point, attention is directed to the attitude angles * and 
a *. These angles have not been employed except in the discussion of the 
steering losses in the material which preceded these discussions (rather, 
combinations of constants were employed); however, it was noted that such 
data could be employed by defining ©<, and /3 terras of a reference value 
and a perturbation as 

ck - cx * y- A 

- /?* / A/2 

This approach was taken in References 1.1 and 1.2. By way of contrast, the 
results do not appear simpler nor do they appear to be more readily mechanized. 
Rather, as was noted in the discussion of the simplified in-plane guidance 
equations, the effect is to distort the similarity in the in-plane and out-of- 
plane control problems. Further, no improvement in the degree-of-optimality 
is achieved, since in both cases the limiting assumption is that the change 
in the attitude during a given burn is small. In way of defense of the 
reference, it is conceded that the approach which he outlines will be subject 
to different numerical problems than the one outlined here (i.e., the problems 
of round-off and loss of numerical significance may be less severe when the 
angles ot* ( /3 ) and AoC. ( A/3) are of different magnitudes). It is 

also conceded that his representation of the coupling between the in-plane 
and out-of-plane dynamics is more general. (The two should agree through the 
third-order terms in the series for cos c< which should be more than adequate 
if the launch time is selected such that the desired plane-of-motion is 
nearly attained) . 

After several values of time-to-go have been predicted in this manner, a 
correction cycle can be superimposed on the result to provide improved pre- 
dictions early in the flight by employing the memory afforded by the preceding 
times. This fact will assure improved accuracy in the generation of the 
steering command and better efficiency in the utilization of the propellant 
available. However, at this time no experimentation has been conducted to 
indicate the best form for this correction cycle. 

Cherry in References 1.3, 1.4 approaches the problem of finding the 
time-to-go in a different manner. His approach is based on the assumptions 
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that the rate of change of the angular momentum per unit mass of the vehicle 
is the moment of the applied force per unit mass (i.e., gravity is assumed to 
be along the instantaneous radius vector) and that the steering angles oC. and 
a are <<1. Thus 


H CO 


ret) x F Ct) 


^ rCt) FCt) cos/3Ct ) 


r Ct) f Ct) 


/3 a Ct)~ 

2. 


Now if an expansion for r = r (t) is assumed, this equation can be integrated 
from 0 to T and the result equated to the desired change in angular momentum. 
The problem with this approach is that any expansion for r = r (t) which 
matches the two ends of the arc will involve the time-to-go (see the dis- 
cussions of average gravity). Thus, it is necessary to utilize a less precise 
representation in order to evaluate T and then to resolve the problem 
utilizing a more precise series. Assume 
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This equation can now be solved iteratively by assuming that ( A v / T 0 r 1 
to obtain a first estimate for T 
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and then by modifying the process to utilize the complete equation 
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This latter equation can be iterated (solved by the method of false position) 
until convergence to a desired degree is obtained. Further, if it is desired, 
the final value of T can be employed to evaluate the constants of the linear 
tangent steering program for /3 and the correction to T resulting from 
steering can be produced by integrating 
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dt 


and resolving for T. However, this form of the correction cycle is con- 
sidered too involved to be practical. Rather, a numerical correction cycle 
is considered preferable if Cherry's method is to be employed and if a 
suitable numerical technique can be devised. However, since more assumptions 
are required to obtain reliable estimates from this formulation, and since 
there could be trouble in the iterative solution process, the former method 
is preferred as an approach to determining the time-to-go. 


2.1.5 Determination of the Time Average Gravity Vector 


As the final step in the derivation of analytic steering equations for 
a single stage vehicle, it is necessary to construct a simplified model of 
the equations of motion. This process in turn is accomplished by considering 
the definition of the time average gravity vector. 


where 



(2.1.54) 
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This process requires that g be expressed as an explicit function of time. 
Consider the case where g will be approximated by an expansion of the follow- 
ing form 
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where A is a vector constant used to assure that g(T)- g(rf) . But, 
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Consider the X component of gravity 


fxf 


~u/Z 


Differentiation yields 



Similarly for the y and z components. The result is 
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where I is a 3 x 3 identity matrix. Thus, since this matrix is known, the ^ 
corresponding vector d/dt (g)/ 0 is known. This fact allows the vector A T 
to be evaluated as 
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Note that the product A 


is well defined for all times 


0 < T £ T, 
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Finally, the time average gravity vector can be computed by substitution 
of equations 2.1.55 into equation 2.1.54 and integrating. 
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2.1.6 Extension to More Than One-Stage Vehicles 




The analysis of both the in-plane and out-of-plane guidance problems 
showed that the steering equations to be mechanized were independent of the 
time dependent nature of the thrust. This fact means that such variations as 
changes in the thrust due to progressive burn of a solid propellant, or dis- 
crete changes due to staging the vehicle do not alter the steering considered 
optimal. Thus, the procedures outlined on the previous pages can be employed 
without change for multi-staged vehicles of differing thrust levels, provided 
that means of estimating the burn time for each of the stages can be formulated. 

To accomplish this expressed objective, it is assumed that the various 
stages have been designed in a near optimum manner, and that their selection 
for a particular mission was based on criteria of optimality which need not 
be considered here. It is also assumed that the propellant loadings for the 
various stages were selected based on rational logic. Under these assumptions, 
the optimal policy for utilizing the propellant is to expend all of the pro- 
pellant in a particular stage (adjusted for loss and ullage) before igniting 
the next. This being the case, the maximum (average) burning times for each 
of the stages can be computed as 


Kyi (2.1.59) 


At this point in the process, the time-to-go is defined as 
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where T n denotes the unknown duration of burn of the last stage to be employed. 
Now, as before, an estimate of the quantity Tg 0 is generated by solving 
equation ( 2 . 1 . 41 ) and the velocity constraint equation (2.1.48) dimultaneously. 
The equations for this solution are: 


and 
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where A\f± (i=j...n-l) denotes the ideal increment in velocity corresponding 
to the maximum burning times for each of the stages below the one of unknown 
burning duration; where j denotes the number of the stage which is presently 
burning; where AV C is the correction or adjustment in the AV obtainable 
from the jth stage to compensate for any ellapsed burning time for the stage; 
and where AV n is the required increment in velocity for the nth stage. 

The procedure employed to solve this equation is as follows : 

1) assume the number of stages n [this step will not be required once 
the first estimate of T n £ T nmax is obtained (unless subsequent 
computations of T n are larger than T njnax ) since the last estimate 
of the number of stages can be employed} . 

2) solve equations (2.1.60) and (2.1.61) for T n in the same fashion 
used in the single-stage formulation (in the single-stage case 

v eff = l^f - V 0 ) + V c ). 

3 ) check T n against the maximum burning time available for the stage. 

If T n is negative, reduce the number of stages by one and repeat 
the process. If T ff is greater than T nmax , increment the number of 
stages by one and repeat the process. If 0 < T n < T njn , the 
number of stages and the time-to-go estimates have been defined 
and the steering constants can be computed. 


The danger in this process is the possibility that the initial estimate 
of g aV g will become much less accurate as the burning times increase, due to 
the fact that the change in the radius was assumed to be reasonably small in 
order to yield a solution for the average gravity. This problem can be 
alleviated by defining an average gravity for each stage. However, such a 
step requires that intermediate radii be established (iteratively) for the 
terminal position for each of the stages. One means of accomplishing this 
objective would be to employ the estimation of g avg presented in equation 
(2.1.58) to define a trajectory [? - ? (t)] ; and segmentation of this tra- 
jectory at times corresponding to the stage burn times to provide a means 
of generating improved estimates of the average gravity for each of the 
stages. (This process should converge quite rapidly, since gravity Is not 
the dominant acceleration.) The term 



can then be replaced by 


-n 






No numerical experimentation has been performed for this case. Conse- 
quently, no estimate of the necessity of the iterative approach is readily 
available. Rather, it is noted that unless the problems are severe, the 
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mechanization difficulties and the smal 1 improvement obtained for most tra- 
jectories of interest would dictate that the approach presented for single 
stage vehicles be employed. 

2.1.7 Extension to Fixed Time Coasts 


No assumption has been made in the previous analysis which would pre- 
clude calling the stage a coast; i.e.. 



0 


To - 7" 

* Const £ 

Thus, if the durations of the coasts are added to the total time in the 
computations of the gravitational loss, and if the velocity increments for 
the corresponding "stages" are equated to zero the result is 

-%-K [Er< .£ r c ,r„] 

^ - AV C+AVn] -/ [0 TJ TTi] 

The probability of error in the estimation of g a is now reasonably 
high. Thus, it is recommended that the iterative logic suggested in the 
extension to multi-stage vehicles be employed or that intermediate radius- 
velocity terminals (r, v will define the coast arc), be established so that 
the total problem can be solved sequentially. This latter alternative is not 
unreasonable since the duration of the coast periods can be rationally selected 
only by performing a study of the effects of these durations on the requirements 
for the boost vehicle. Thus, if a fixed duration coast is to be commanded 
(based on simulations of the trajectories, etc.), it is reasonable to assume 
that the "terminals" can also be supplied. If so, each phase of the problem 
(K consecutive burns) can be considered as a separate problem. 

The alternative to requiring intermediate terminals (thus losing flexi- 
bility in the targeting for the vehicle) is a more complete simulations 
capability for the guidance computer, for the crew and/or for the ground based 
tracking station. 

2.1.8 Results of a Typical Simulation 

The material of Sections 2.1.2 through 2.1.5 has been mechanized for 
numerical simulation to demonstrate the nature of the solution and the type 
of control derived and to define the general nature of the mechanization re- 
quirements. This simulation is illustrated in Figures 2.1.1, 2.1.2, 2.1.3, 
and 2.1.4. Figure 2.1.1 shows an overall logic for the guidance function and 
demonstrates the effects of a time delay in generating new guidance commands 
on the system mechanization. Figure 2.1.2 presents the first step in the 
process, that of the determination of time-to-go as a function of the error 
signals in position and velocity Once the first estimate of time-to-go is 
made, estimates of the gravity vector and the steering constants generated 


35 



utilizing the logic of Sections 2.1.5 and 2.1.4 as illustrated in Figures 2.1.3 
and 2.1.4. At this point, the gravitational and steering losses are estimated 
and a correction to the estimated time-to-go (two corrections are made the 
first pass through the process). At this point, a corrected time-to-go can be 
predicted from previous memory, if available (this logic has not as yet been 
checked). The revised time estimate is then utilized to predict the guidance 
steering constants. 

One point of particular merit should be noted here in regards to the comp- 
utation of the steering constants . The possibility exists that the estimated 
time-to-go will be so poor that the steering constants C and K]_ (representing 
the sines of the angles a and g at the epoch t 0 ) will be larger than one. If 
this situation exists for either constant, the corresponding steering angle(s) 
should be set at 90° (depending on the sign of C and Ki) until the problem 
becomes better behaved. 

Once the steering angles are defined, the motion of the vehicle is simu- 
lated for a small interval of time (the interval for the simulation is m times 
the step size of integration, where m is an input quantity intended to show 
the effects of appreciable computation times, i.e. , time delays, on the 
resultant trajectory). The process is then repeated until the time-to-go 
reaches negligible proportions of until the magnitude of the velocity matches 
that which is required. 

The results of two particular simulations are shown in Figures 2.1.5 and 
2.1.6. These figures were generated for the problems: 
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respectively. The figures show the manner in which the position relative 
to the desired injection point goes to zero and the nature of the approxima- 
tions for the time-to-go and the steering constants. The figures also show 
the nature of the gravitation approximation by defining the constants of 
proportionality necessary to construct the "average" gravity vector from the 
instantaneous gravity vector and the terminal gravity vector. Note that for 
the major portion of the trajectory, the components of the average gravity 
in the g D direction are nearly constant. (The sharp increase in this factor 
late in the flights is the result of indeterminancy in the solution for the 
proportionality constants as the elapsed time between t Q and t^ goes to zero, 
i.e., go -*> gf. No problem was encountered in the definition of the gravity 
vector itself). This fact implies that studies of the gravitational model 
employed in the analysis might reveal an extremely simple yet highly accurate 
empirical relationship which could be used to replace the present computation 
of the time average gravity vector. 

Figure 2.1.7 demonstrates a measure of the off-optimum nature of the 
solution. Presented is a plot of the first estimate of the time-to-go for a 
series of problems identical to the samples which were discussed but with 
different components of velocity in the R direction at the time guidance is 
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initiated. Also presented is a plot of the "true" time-to-go for the same 
runs. Note as the amount of steering increases (zero steering occurs for 
r 0 ^ 9000 fps), the initial estimates become increasingly poor due to the 

assumptions made in developing the estimate. In all cases, however, the first 
estimate was better than that obtained without consideration of the steering 
losses. For the purposes of contrast, this latter estimate and the corre- 
sponding true time-to-go is also presented. Of particular interest is the 
fact that the initial guesses without the steering losses are valid (under an 
assumed accuracy constraint, e.g. 10$) over a smaller range and that more fuel 
( T true^ ls acquired. 

These results by no means exhaust the evaluation of the adaptive (itera- 
tive) guidance scheme. They are intended to show the nature of the solution, 
and to demonstrate possible revisions which will make the process more accu- 
rate, more efficient, and/or more flexible. Thus, continued effort in this 
area will be followed with great interest. However, the simulation does 
serve to indicate the approximate nature of the mechanization requirements 
(no attempt has been made to optimize the FORTRAN coding or core storage as 
would be the case in a true self-contained system) . The resultant require- 
ments are; 

Time-to-go 1000 octal (including steering losses) 

Gravity 400 octal 

Guide 500 octal 

1900 octal 

In addition, the following general purpose subroutines were mechanized: 
square root, natural logarithm, sine, cosine and dot product. The exact 
length of these routines is not known. 
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Figure 2.1.1 

Iterative Guidance Mechanization 












Figure 3.1.2 
Time-To-Go Estimation 
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Figure 2.1.3 

Average Gravity Computation 
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Figure 2 . 1 . 1 + 

Steering Constant Computation 
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Figure 2.1.5a 

Terminal Motion for Iterative Guidance 
(Sample l) 
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Figure 2.1. 5d 

Out-of-Plane and In-Plane Thrust Attitude 
Angles for Iterative Guidance 

(Sample l) 
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Figure 2.1.6a 

Terminal Motion for Iterative Guidances 
(Sample 2) 
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Figure 2.1.6b 

Predicted Time of Flight for Iterative Guidance 
(Sample 2) 
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Figure 2.1. 6d 

Out-of-Plane and In-Plane Thrust 
Attitude Angles for Iterative Guidance 


(Sample 2) 
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2.1.9 Extension to Constrain the Longitudinal Component of Position at 
Arrival 


The discussions presented in the previous sections have been directed 
toward the formulation of the boost guidance problem in which only the radial 
and lateral rates and the terminal radius and lateral position were constrained. 
The present discussion is intended to introduce an approximate means of con- 
straining the longitudinal displacement and rates. 

Consider the equation for longitudinal motion 
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where cos is assumed to obey equation ( 2 . 1 . 30 ) 

sino^ is assumed to obey equation (2.1.35) 
and consider the approximation of cos /§ by employing the identity 
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This equation can be utilized to generate the solution for S(T) and S(T); 
however, the solution will be non-linear in the constants Ki and K2- 
Further, there will be a set of four (4) boundary conditions which must be 
matched by selecting these two constants. Thus, in general, the equations 
will not yield a unique solution even under the assumption that they can be 
easily solved. This fact has lead to the partial reformulation of the 
problem around the lines suggested in equation (2.1.29). That is, additional 
constants will be introduced which will be selected so as to satisfy the 
terminal constraint on longitudinal motion. 

Consider the equation (2.1.33) and (2.1.34) which can be expressed in 


the form 
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These linear approximations are required since the solution to be performed 
on the guidance computer must be linear in the steering constants to avoid 
iterative computations. Note that while terms involving t^ could be easily 
added to the representations of both the sine and the cosine, there would be 
an insufficient number of boundary conditions to provide the constants with- 
out employing the dependence of these functions (at prescribed epochs) in the 
form 


i 

sin 


/3 
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a 

Co* 


/3 


/ 


However, since this identity has not been employed, it is highly pro- 
bable that the resultant equations will be valid only for relatively small 
ranges of initial and terminal states. No attempt has been made to determine 
if this limitation exist or the severity of its effect on the steering program. 
One saving grace exists though since the results of the sample problems in-, 
dicate that the steering angle /2 is quadratic (roughly) in time. Thus, 
is roughly linear and 

/ 3 ) 
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Since this result is the same as that obtained in the previous series of 
approximations, it is assumed that the result is sufficiently accurate to 
allow for an approximate solution to the equations of motion. 
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Now, under the assumptions that these approximations are valid for some 
trajectories, the steering constants K^, K^, Kc and K/ can be evaluated by 
employing the solution presented in equations (2.1.37) and (2.1.38) by simply 
changing notation. 

The result is 
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where 


Ar - r-CD- r(0)-f Avf * £ 7 
zir = r{r; -r(0)-f n J' * 0) T 
AS = S(T)-S(0)- ? ; 7 *ST 
ziS « S(T)-S(0)'^ v/ -S ^-S(0)T 

A u > A JJL , A XI f A tl _ defined in equation (2.1.39) 

Finally, since the coefficient matrix relating the steering constants is 
partitional as it is, the solutions for the constants is 



and the "optimum" steering angle is given by 


Taw /S 


<, - k v t 

Ks ~ *1* 


Note that because of the previous approximations, the solutions for the sine 
and cosine of (3 are independent and that the constants and are 
identically and respectively. This fact allows the approximate con- 
straint on the longitudinal motion to be added or deleted from the guidance 
solution as desired without modification to the formulation by simply bypass- 
ing a portion of the logic. Care must, however, be exercised to assure that 
the assumptions implicit in the development are not violated. 
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2.2 EXPLICIT GUIDANCE EMPLOYING GUIDANCE POLYNOMIALS 

2.2.1 Introduction 

The development of rocket vehicles that are capable of Injecting multi- 
ton payloads into orbit has established the need for a guidance system that 
differs from those that were developed for ballistic missiles. Further, 
the high cost of the vehicles, along with their utility in the Apollo program 
to carry astronauts makes it necessary to have a guidance system with a 
high degree of reliability. These objectives can be attained by increasing 
the flexibility of the guidance system to make it compatible with the charac- 
teristics of these vehicles. For example, the enormous thrusts that are 
generated by the early stages are normally achieved by clustering engines. 
Therefore, if one of the engines fails a discrete variation in the thrust 
would result at an unpredictable point on the trajectory. This failure would 
not, however, effect the total energy available to complete the flight but 
would result in a lower rate at which energy is available; this fact would 
in turn alter the shape of the trajectory which can be flown. However, if 
the guidance system has been designed to accommodate this type of failure, 
it is still possible to complete the specified mission. 

Further, the large vehicles are expected to be used for a variety of 
different missions. Thus, a guidance system with sufficient flexibility 
to be used for each of these missions without major redesign or modification 
is highly desirable due both to considerations of cost and reliability. 

And finally, the complexity of these vehicles are such that the probability 
of their achieving lift-off at a specific instant is small. Thus, the gui- 
dance system must be designed to compensate for these variations. 

The objective of this section is to explore one approach to the problem 
of providing a guidance logic adequate for all of these requirements. To 
this end, the following paragraphs have been prepared. 

2.2.2 Preliminary Considerations 

In the discussion to follow, it will be assumed that a navigation system 
is available that indicates the vehicle position, velocity, acceleration, and 
attitude in a continuous mode. For example, such a system could be either 
inertial or radio navigation. If an inertial navigation system is used, then 
accelerometers are mounted on a gyro stabilized platform and the position 
and velocity of the vehicle are computed from the accelerations that are 
measured using these instruments. If the radio guidance system is used, 
one or more radar units will observe the distance of and direction to the 
vehicle, and this information will then be used to compute the position 
velocity and acceleration. 

The guidance process is defined as the plan by which the navigation 
information is used to control the flight of the vehicle. For a rocket 
vehicle this process will take the form of two sets of equations. The first 
set is defined as the steering equations and is used to compute the direction 
the thrust vector should have for the vehicle to achieve the desired flight 
path. The second set of equations is used to compute the engine throttle 
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setting. In general, this seeond set includes (l) the time of thrust initia- 
tion, (2) instantaneous thrust magnitude, ( 3 ) time of thrust termination. 
However, the large liquid vehicles presently being built have fixed thrust 
engines so that the instantaneous thrust magnitude can be controlled only 
in a small region about the nominal. As a consequence, only the times of 
thrust initiation and termination is generally computed. 

A number of different guidance processes have been developed. Among 
this set, are those for which the steering equations have been intuitively 
selected. For example, one of these schemes orients the engines so that the 
thrust acceleration is in the same direction as the velocity-to-be-gained. 

(The difference in the velocity vector required to achieve the specified 
terminal position at some future time in free flight and the instantaneous 
velocity vector). In this system the navigation system is used to monitor 
the vehicle flight, and the thrust is terminated when the velocity-to-be- 
gained has been driven sufficiently close to zero. The time of thrust 
initiation is selected to give some desired end condition, with the aid of a 
simulation of the process. It is probably possible to freeze a process of this 
nature and engineer extensions that will handle particular missions. However, 
the previously mentioned objections still persist, the (l) the knowledge 
of the applicability of a guidance process is limited to one mission; (2) 
the degree to which the process reaches its theoretical optimum is not known. 

For reasons such as those mentioned, the decision was made to start 
afresh with new concepts that were independent of vehicle and mission 
configuration. The objective of this new approach was to develop a guidance 
process which would adapt to the particular vehicle and mission being flown 
(thus, the designation path- adapt ive ) . In this approach, the same basic 
guidance plan will be applied to all missions and configurations, with the 
steering and throttle equation being selected to minimize a particular loss 
function. 

2.2.3 Formulation of the Guidance Equations 

In order to gain some insight into the areas over which an optimization 
of the guidance mode may be made, consider Figure 2.2.1 
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In this figure, the flight space Is assumed to he two-dimensional, and the 
objective is to achieve the mission trajectory. The mission will be repre- 
sented mathematically by the set of equations 

fi * ) = 0 1,2,3 (2.2.1) 

where the solution to the system of equations in x, y, x, y, and t Is, 
of course, not unique, since any point on the curve represents a solution. 
Therefore, rather than choose one particular solution as a standard injection 
point, the mission is left specified in a functional form. This degree of 
generality assures that guidance may be generated at any given point in the 
flight, /l o , which is which is "best" in the sense of some scalar com- 
parison function and which will satisfy the mission. 

If desired, an arbitrary form of thrust direction history could be 
assumed, and the corresponding optimization for the best solution of the sys- 
tem of mission equations (2.2.1) made. It seems wisest, however, to seek 
that form which is optimum over all possible functional forms and to use the 
one function of that family which also provides the optimum solution to the 
mission equations. That is, among all steering functions which result in 
flight paths that have end points that satisfy the mission conditions, the 
one that minimizes some specified quantity is sought. This type of problem 
is treated in the Calculus of Variations, where theory exists that is useful 
for singling out the one particular desired function. 

2.2.h Sample Steering Function Solution 


In order to obtain a guide to the approach and methods necessary to 
implement the adaptive guidance mode, using the concepts of the Calculus of 
Variations, it is helpful to consider a simplified problem; thus drawing 
attention to those steps that are characteristic of the methods applied to 
actual flight problems. 

The particular problem considered was selected for two reasons. First, 
a closed form solution can be given for it. Second, the simplification of 
the problem allows a closer exposition of the salient features of the approach 
taken. 


Consider the flight of a vehicle in an inverse square gravitational 
field that is propelled by a rocket having a constant thrust and mass flow 
rate. This physical situation can be represented by the following system 
of equations. 





( 2 . 2 . 2 ) 


si = A t - U* si 3 y? 


(2.2.3) 
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where 


A = • /i =■ A j. 0 ~ ~ X 

St 

(2.2.4) 

K z 

y = -4 ~A r C*o9- — if 

(2.2.5) 

Assume the initial conditions 


, \ 
’ -- ’ -<■<■ • 

(2.2.6) 

to be ’known and that the mission criteria expressed by equations (2.2.1) are 
to be achieved with a minimum thrust time ( - t 0 ) . The steering function 

which does this can be determined via the calculus of variations to satisfy 
the following equations 

*y] 

(2.2.T) 

~ 7 Zf + ( £Ly* 

(2.2.8) 

0 - A, c&o 9 - \ z 9 

(2.2.9) 


(where the initial conditions for the Lagrange multipliers for each component 
of the terminal state which is unconstrained may be selected arbitrarily and 
where the remaining multipliers must be determined by trial and error) 
together with equations (2.2.2), (2.2.4), (2.2.5), and the end conditions 
expressed by equations (2.2.1) and (2.2.6). 

This set of equations constitutes a two-point boundary value problem. 

One form of the solution would be the set of functions QU), Alt) > A It) 
that meet the stated requirements, based on the absence of any disturbances 
in the interval from if to i a . If a solution is obtained based on a 

state which is taken later on along the mission trajectory, the same set of 
functions yi(t) , Act) and Bit) would be obtained. 

However, as soon as a disturbance occurs, a new set of solutions will 
result. Thus, the solutions Alt) , Alt) , and Git) are in turn functions 
of the initial conditions and could be written explicitly as functions of 
these parameters. In particular, the function Bit) evaluated at 
t - i a is 

= (2.2.10) 

*“ 7ft O J 
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This equation expresses the solution of the desired thrust direction in terms 
of the measured state at the initial, epoch. The remaining end conditions, 
while they are not given explicitly, can be expressed similarly. 

The resulting set of equations for the unknown end values in terms of 
those that are known may be thought of as a second form of the solution to 
the two-point boundary value problem. The most important of these, however, 
are equation ( 2 . 2 . 9 ) and a similar expression for if . Thus equation 
(2.2.9) is the steering function and the expression for if would be the 
thrust termination function. 

The success of path-adaptive guidance is determined largely by the qua- 
lity of the equations that are mechanized in the guidance computer. This, 
in turn, depends on the nature of the function necessary to accomplish the 
desired results and upon the manner in which it is prepared for the computer. 
However, guidance functions that represent the optimum exactly (to measurable 
accuracy) can be expected to require more computer complexity and weight 
than simplified approximations to these functions that produce approximate 
opt imums . Thus, trade-off studies are required to determine the allowable 
degree of approximation for each of the applications. 

2.2.5 Analyti ca l Approach to Guidance Function Representation 

The form of the solution of the two-point boundary value problem dis- 
cussed previously (equations (2.2.l) to (2.2.9))* may be attempted analyti- 
cally. Concept ually , this solution could proceed as in the following 
simplified example. (Analytic solutions to more complicated problems are 
discussed in Reference 2.1 and 2.2) 

Consider the problem of steering to a point in the ( V ) space from 
arbitrary initial conditions 


U n 


0 •> 


V 0 


( 2 . 2 . 11 ) 


to the fulfillment of the mission criteria. 


If the equations of motion are 



j/f - // 0 

(2.2.12) 


Vf ~ V ~ 0 

(2.2.13) 


= S/AJ OC 

(2.2.14) 

O 

V 

r eos oc 

( 2 . 2 . 15 ) 


and if the objective of the guidance system is to minimize ±p - i D , 

application of the calculus of variation provides the equations that define 
oc as 
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A = 0 

(2.2.16) 

-4 

(2.2.I7) 

£OS e*c - jgg S/Me* = O 

(2.2.18) 


These equations have the solution 


(2.2.19) 

Thus, upon substitution of this solution into the equations of motion and 
integrating from t 0 to tg 

-^U c = (tf ~ * 0 ) S/A/ c < 0 (2.2.20) 


- Vo - (*-f - ■to') £OSoC 0 ( 2 . 2 . 21 ) 

This system is then solved for //p and V„ substitution made into equations 
(2.2.12) and (2.2.13). f / 


AJf. -^U a = C*4 0 ) S/A/oZc 

(2.2.22) 

V^-Vo ^(±4 -to) eos oc 0 

(2.2.23) 

The simultaneous solution provides the steering function 

XLS - 

<y- 0 - CUUI toy, t — - 

(2.2.2k) 

Vf - Vo 

and the cutoff function 

- i *Vc u f -ju.)* + (v f -Vo) z 

(2.2.25) 

Now consider the more realistic problem discussed in 

section (2.1) in 


its most complete or in its simplified form. For this problem, the equations 
of motion are more involved and the variables are coupled with the result 
that no simple solution for the Lagrange multipliers or the corresponding 
steering angle can be obtained. Thus, the analytic solution (or simple 
iterative solution) required for an onboard guidance computer cannot be 
realized. This fact has lead to an empirical approach to the problem; this 
approach will be presented in the following section. 
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2.2.6 Guidance Function Generation Using An Empirical Approach 

For this approach, the two-point "boundary value problem represented by 
equations (2.2.1) through (2,2.8) is solved numerically for a particular assumed 
set of Initial conditions, equation (2.2.6). This process is repeated for 
a large variety of values that lie within a region that contains all of the 
disturbances that the vehicle is designed to withstand. Thus, a tabular 
representation of the guidance functions is constructed which may then be 
approximated by a polynomial. This task has thus become a curve fitting 
problem for a function of several variables. Wo general theory presently 
exists to accomplish this objective; however, some relatively good results 
have been obtained using intuitive procedures. 

For evaluation of this approach, a test problem was postulated. The 
task was to derive a steering function for the vacuum flight of the second 
stage of a vehicle. Thus, the initial point is the cutoff point of the 
first stage. Wo cutoff function was employed, rather it was assumed that 
the fuel would be burned to depletion. (The application of this technique 
to other tasks is discussed in References 2.4 and 2.5*) 

The procedure used to build the table for the steering function was to 
isolate a family of optimum trajectories originating from the area of the 
first stage cutoff conditions and satisfying the terminal end conditions. 

The family consisted of 126 trajectories for which variation had been made 
in the upper stage thrust level. It should be noted that each point on any 
one of the trajectories was in fact another initial conditions, thus 
providing another value for the table. For each trajectory of the family, 
points were read every five seconds to provide a time history of the steering 
function. 


*o[js jfe jdLr* ] ( 2 . 2 . 26 ) 

Since this mission is independent of time, does not occur in the steering 
function. 

The curve fitting procedures chosen to represent (2.2.28) was the method 
of least squares, since it was felt that this process would lead to a polyno- 
mial form that is especially convenient for an on-board computer. Further, 
the polynomial that was used to approximate equation ( 2 . 2 . 26 ) was of the form 


& j = W ‘J T W *j f a 3^J “• W ij 
where W-j are generalized product functions of the type 


( 2 . 2 . 27 ) 




h 



( 2 . 2 . 28 ) 


and where the indices range over the set that includes all postulated powers 
for the various factors in the series (h, k, p, q, r and s; generally these 
powers must he assumed, the guidance function determined, flight simulated, 
and the exponents empirically optimized by trial and error). 
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The choice of the particular powers for the variables to be contained 
in the set of indices will depend on the function being simulated. A method 
for selecting these indices is not yet available, however, for any particular 
choice it is desired to find the best set of coefficients, olj . Thus, at 
this point it is assumed that the polynomial has been empirically defined and 
that the problem is now to define the best set of coefficients in the sense 
that the sum of the squares of the differences between the values provided 
by the polynomial and optimum values of © are a minimum; i.e., 


L = £ (2.2.29) 

is a minimum (where, © , is the optimum value of the thrust attitude). But 

the computational algorithm for defining the a's is 


» \mm t \ mq‘ 


(2.2.30) 



a, 


e, 


a z 

Jxi 

£ 

a 

2 <2, 

9 ‘ 2 





/ 

9 - 




j 


X 

w lZ . . . 

w tj 

t*j 

w z\ 

w zz . . . 


M 5 

X 

Wst ■ • • 



X. 


W *J- 


Note that the index j ranges over all of the tabulated points and that the 
number of tabulated points exceeds the number of coefficients, i.e., 

(This latter requirement assures that a solution can be generated) 


By varying both the selection of the tabulated points and the specific 
polynomial form used in the fitting, the residuals can be modified so that 
the loss function of equation (2.2.29) is as small as possible. The final 
criteria of goodness of fit, however, must be a weighing of the number and 
complexity of terms in the polynomial against the degree of optimization and 
accuracy of meeting mission conditions achieved by the polynomial. For 
fitting the 12,600 point tabulation, various approximating polynomials up 
through third order were investigated, (Reference 2.3) and it was found that 
a third order polynomial of eighty-four terms proved the best compromise 
representation. This conclusion is predicated on the observation that the 
errors in the control deflections relative to the true optimum policy were 
of the order of .15 degrees while assuring that the storage limitations of 
the on-board computer were not exceeded. 
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In conclusion it is noted that while the entire discussion of this sec- 
tion of the monograph has been slanted toward the two dimensional problem, 
the results are valid for the three dimensional problem as well. The solution 
will in all cases be more complex due to the fact that additional terms in 
the equations will be required (unless the in-plane and out-of -plane motions 
are separated as discussed in section 2.1). However, the concepts are 
unaltered. 
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2.3 PERTURBATION GUIDANCE 


This section of the monograph will deal with a class of guidance 
schemes which will be referred to as perturbation guidance. This name has 
been selected because each of the schemes in some way relies on the closeness 
of the true flight to precomputed information that has been computed for a 
specific reference flight. In other words, a significant amount of computa- 
tion is done on the ground prior to the flight for a particular mission in 
order to minimize the amount of on-board computation. Since the preflight 
computation is performed with a particular mission in mind and with a nearly 
exact model for all forces acting on the system, it is assumed that the 
preflight results will be in the "neighborhood" of the actual flight results. 
It is this "closeness" of the preflight computation that is the basis of any 
perturbation guidance scheme and which affords accuracies comparable to those 
obtained if a large and very fast computer was used for real time computation 
during the actual flight. As long as the actual flight does not differ signi- 
ficantly from this "nominal" or reference trajectory, the precomputed infor- 
mation is sufficiently accurate to allow the guidance scheme , which is 
designed by approximating some quantity by a Taylor series, to correct for 
measured deviations in position and velocity. The original assumption of 
closeness to a nominal trajectory permits the Taylor series to be truncated 
after first order terms in most cases without introducing intolerable errors. 
This approximation then permits the use of the more powerful techniques of 
linear analysis in many cases . 

Although the techniques presented herein have many similarities, certain 
peculiar features have been singled out in order to distinguish the schemes. 

It should be noted, however, that any attempted organization based on these 
peculiarities is artificial since similar mathematical techniques are used in 
conjunction with al 1 of the guidance schemes. In this light, the organiza- 
tion of the section will now be presented. These discussions are divided 
into two major parts. The first part, 2.3.1, deals with guidance schemes 
that use the required velocity concept (defined in text). The second part, 
2.3*2, presents a class of guidance schemes that reduce the nonlinear equa- 
tions to linear perturbation equations . Following the derivation of the 
equations for each of these guidance schemes, a control or steering section 
is presented. Since it is usually desired to optimize the control in some 
sense, the determination of the optimum control policy involves the use of 
calculus of variations, Pontryagin's Maximum Principle, and/or dynamic 
programming. A sample application of each of these approaches is included. 
Section 2.3*1*2.3*1, Optimum Steering for C* Guidance, is an application of 
Pontryagin's Principle. Sections 2. 3-2. 2.1 and 2. 3*2. 2. 2 are respectively 
applications of calculus of variations and dynamic programming. It should 
be noted that a detailed discussion of the variational techniques used in 
these sections is beyond the scope of the monograph. However, all of the 
required information may be found in other monographs of the series 
(references 3*22, 3*23 and 3*24). 
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Before continuing with a detailed analysis of the various Perturbation 
Guidance schemes, it is worthy to note some of the features of each so that 
the advantages and disadvantages of each can be kept in mind as they are 
presented. Figure 2.3.1 is a table that compares and contrasts the three 
main types of Perturbation Guidance schemes. Delta Guidance, C* (or Q) 

Guidance, and Linearized Perturbation Guidance. These three schemes are 
compared on the basis of performance, optimization, errors, application and 
mechanization . 

2.3-1 Required Velocity Approach 

The guidance schemes in this section share the use of the required velo- 
city concept. (The required velocity is defined as the velocity that is 
needed by a vehicle in order for it to reach a specified position at some 
specified future time under the assumption of free flight from its present 
position.) The two schemes presented in this section are two different 
mechanizations of the equations for the required velocity. The idea behind 
both of these methods is to provide a simple means of updating the required 
velocity as the vehicle progresses along the powered flight. By measurement 
of acceleration of the vehicle the first concept. Delta Guidance, expands 
the required velocity in a Taylor series about some nominal burn-out point . 

A good approximation of the required velocity for a point other than the 
nominal burn-out point can then be found by substituting the position coor- 
dinates of the vehicle into the Taylor expansion. The position of the 
vehicle at subsequent times is then determined by integrations of the total 
vehicle acceleration. This scheme provides a continuous knowledge of the 
required velocity for purposes of steering and is capable for compensating 
for errors at earlier epochs. 

The second guidance scheme that used the required velocity concept is C* 
(sometimes called Q) Guidance. This scheme is a means of continuously 
updating the velocity-to-be-gained (the difference between the required 
velocity and the current velocity.) Measurement of the thrust acceleration 
provides the information that is necessary to continuously update the 
velocity-to-be-gained such that the vehicle can be steered properly. 

Acknowledgement is given to C. W. Sarture (reference 3-3) whose material 
on Delta Guidance was of significant assistance in the preparation of this 
section of the monograph. 

2. 3. 1.1 Delta Guidance 

2. 3. 1.1.1 General Discussion . The object of Delta Guidance is to provide 
a reasonably accurate value of the required velocity throughout all phases of 
a powered flight. However, since the expressions for the required velocity 
are quite complex and non-linear, an ''exact" knowledge of the required velocity 
as the vehicle progresses in its powered flight would require a tremendously 
large and fast computer performing real time computation. To remove this 
problem, it has been found that a Taylor series expansion as a function of 
position and time about the nominal burnout point provides a sufficiently 
accurate value of the required velocity for any position or time along the 
nominal trajectory and at the same time produces a significant reduction in 
on-board computation. 
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Figure 2.3.1 A Comparison of Perturbation Guidance Schemes 
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It should be noted that while such an expansion has advantages from the 
point of view of real time computation, it has corresponding disadvantages in 
that there is no guarantee that the vehicle will adhere to a nominal trajec- 
tory, especially at the beginning of the guided flight where the vehicle is 
its farthest from the burnout point (Taylor expansion point). It should also 
be noted that since most computation is done prior to the actual flight, 

(the results of these computations are the Taylor series coefficients), any 
mission change requires complete pre-flight reprogramming including, at a 
minimum, a new set of Taylor coefficients. 

In general, the expressions for the required velocity are quite complex 
and cannot be written explicitly so that differentiation and expansion in the 
Taylor Series is possible. A grossly simplified problem employing a flat 
earth can be analyzed in order to demonstrate the theory involved. Following 
the flat earth analysis, the more practical numerical techniques that are 
utilized in the preflight computations will be discussed. 

2. 3. 1.1. 2 Derivation of Equations 

2.3-1-1-2.1 General . The required velocity expression in the most 
general case will be a function of position and time, i.e., = F(x,y,z,t). 

Each component of Vp can be expressed as a Taylor Series of a~function of 
four variables as fallows; 




21 


Ax — + Au — + A a — +At — 
dx ? dy <t d/^ H 


+ • ‘ * ■ 


(2.3.1) 


The x component of Vj, then becomes 


= \ * Z' lz ~ K) * f > ) + u lt ' 1 ’ ) 

,2 3 * . ,2 3 2 . .2 )■ 




r ' '■ V 


d*F . .. . d z F . . . . ? Z F 


67 



d*F d* F d z F 

d*F 


d*F 


d z F 


ZF ( t -* )( ** u w ( r^rF\ » ( r*x*-u 


( 2 . 3 . 2 ) 


T d g F 
+ dt?% 


d z F d z F 

U - 4 )(z-x.)t ^ l i-t. )(?-?, )f dtii U ' i - )i i > \ 


* ... 


where V^ x is the nominal burnout velocity in the x direction. Similar expres- 
sions result for v Ry and v R . As mentioned above, the partial derivatives 
are usually evaluated at tne nominal burnout point (x , y , z Q , t Q ) so that 
the accuracy in the representation improves as the desired terminal state is 
approached. These partials are retained as constants throughout the thrust 
phase. The required velocity then is a function of present position and 
present time. 

2. 3- 1.1. 2. 2 Simple Flat Earth Example . An over-simplified flat earth 

problem employing a uniform gravitational field will now be analyzed so that 
the previous theory can be interpreted clearly. Consider a short range 
ballistic vehicle whose target coordinates are designated by (xip, Vp, tp), 
where tip is the desired time of impact. The free flight rectilinear equa- 
tions of motion for this problem are : 


X r -X,tv x U-0 (2.3.3) 

fr V. (2 ‘ 3 - 4) 

More specifically, the equations for the required velocity for free flight 
target impact at the designated time from any burnout point (x^, y^, t^) are 

X T ~ X. 

v e = 

* ^ ~ t b 

U — jS T ~ <^ b X % ( ~ ) 

*t tir-O 

The partial derivatives for these expressions can now be formed in a straight- 
forward manner by treating the burnout coordinates as variables. 


(2.3.5) 

( 2 . 3 - 6 ) 
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If the nominal "burnout and target parameters are assumed to be 


= 85,000 feet y^ = 123,000 feet 

Xrp = 1,212,000 feet = 0 feet 

then the partial derivatives become 


t-^ = 70 sec. 

t T — 392 sec. 


dv p 

— - -3.10 X/O 
dx L 



dX, 


O 



9V** 


10. 9 



» - 3.1x10 ' 3 



-17.3 


The only constants remaining to be determined are thus the nominal required 
velocities at the nominal burnout point. These can easily be determined 
from equations 2.3-5 and 2.3-6 by employing the nominal burnout and target 
parameters . The result is 


= SSOO /*-/-**■ 

- 4800 

io ' 

Hence, the Delta Guidance equations for the required velocity in this parti- 
cular case become 

\ - 3500 -(3.1 x ZO~ 3 )( x- 85,000) +(/0.9)U-70) (2-3-7) 

V^-4800 -(3J X/0' 3 Xy-J23,700)-(l7. 3)(l-70) (2.3-8) 


The required velocity is used by the steering law in order to insure that the 
thrust is in the correct, direction in order to satisfy the desired terminal 
condition. The following sketch illustrates the guidance loop for Delta Guidance. 
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The discussion of steering policies for this and other schemes is presented 
in section 2. 3. 1«3* 

2, 3- 1.1. 2. 3 Numerical Techniques . In a more realistic problem, the 
equations for the required velocity cannot be found in a form which will 
allow the partial derivatives to be determined analytically. However, the 
numerical values for the partial derivatives at the burnout point can be 
determined by simulation techniques. The name given to the most frequently 
used technique is "targeting", a method of determining the best values for 
the partial derivatives in a least squares sense by equating the Taylor 
Series for the required velocity at some perturbed burnout point to the 
corresponding value of the required velocity for the same perturbed point 
as calculated from the best available equations. If the number of perturbed 
points investigated is equal to or greater than the number of unknown 
constants in the Taylor Series, then a "best" estimate can be found for the 
partial derivatives evaluated at the point of interest (burnout point). 

The required velocity for a point slightly perturbed from the nominal 
burnout point can be mathematically stated as 


% ^ % (x*f*2**) 

(2.3.9) 


(2.3.10) 


(2.3.11) 


Thus, the Taylor series (retaining only linear terms) for the required 
velocity can be expressed as 


Vo = v. 


K 


XX 


(x-x 0 )*K X y{y- 






(2.3.12) 
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where = the required velocity of the nominal burnout point. 


X = 


XX 

dx 

- 


x t 

d v 

!C V = 



dV g \ 

^Xi ~ 

»A 


hem. 


and x Q , y G , z q , t D = nominal burnout conditions. Similar Taylor series can 
be written for Vp and V-^ . It should be noted that there is no limitation 
to the process which would preclude the inclusion of higher order terms in 
the series. Indeed, these terms can be added simply once the coefficients 
are determined on the ground. 

Now, since the perturbed point is assumed to be specified (neglecting 
errors in the estimate resulting from B4U errors, etc.), the only unknowns 
in equation 2.3.12 are the partial derivatives. The problem now becomes one 
of finding the best numbers to use for the partial derivatives such that the 
values chosen form a "best" fit in the least squares sense. 

The classical method that is used on this type of problem involves the 
minimization of the mean square error between the values predicted by the 
Taylor Series and those calculated. Thus, if the error is defined as 


^ x ' y t % i ^ ^xx 

where i represents the i th perturbed point and if N perturbed points are 
investigated (N = the number of unknowns), then the mean squared value of 
the error of all N points is 


e z = e* f c z + 

x Xt XZ 


XAJ 


(2.3.10 


The best least squares choice for the unknown constants can be found by 
setting the partial derivatives with respect to each of the unknown para- 
meters equal to zero. The result is a set of simultaneous equations in 
which the number of equations is exactly equal to the number of unknown 
constants. Expressed mathematically. 
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x " dKxx 

P^xa/ __ q (2.3»15b) 
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Now, since the partial derivatives in these expressions 


are 




XL 


2K 




XX 


!!s = -<*-*) 
ffa , -(,.-|.) 




*7 


?£s - -«,-<;.) 

»*« 


(2.3.15c) 

(2.3.15d) 


(2. 3.1 6a ) 

(2.3.16b) 

(2.3.16c) 

(2.3-l6d) 

(2.3.17a) 

(2.3.17b) 

(2.3.17c) 

(2.3*l7d) 


equations (7.3.16) are seen to be a 
unknowns K^, K xyj *bcz * ^xt * 


set of simultaneous equations 


in the 
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The least squares solution to equations (2.3-16) is 

K x = (8 T WB)~'8 T W(AV- ) 


where 


*X = 


K 
^ x 


XX 

f 


k. 


Zi 


^XTl 






/6 ~ fo £, - i a 


t„-t J 


W = weighting matrix which can be used if desired 



= vector of errors in series representation 

Vp (x,y,z,t) - VL . 

K x *x 


with an identical form of solution for the vector constants Ky and K z - 


2. 3. 1.2 G~ :; - Guidance 


2.3 .1.2.1 General Description 

The C* (sometimes called Q) Guidance scheme is another method that uses 
the required velocity concept discussed in the previous section on Delta 
Guidance. More specifically, this scheme uses a. differential equation in 
terms of the velocity-to-be-gained as a means of updating target information. 
The velocity-to-be-gained is defined as the difference between the present 
true velocity and the present required velocity, i.e.. 




(2.3.18) 


The continuous knowledge of the present velocity-to-be-gained then provides 
the information required for steering. The differential equation developed 
in the following analysis is a second order differential equation in Vg with 
the measured thrust as the forcing function. Hence, the velocity-to-be- 
gained is known as a function of time providing the C* guidance equation is 
forced properly. The appropriate steering for C* guidance will be presented 
in section 2.3 -1.3 • 
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2. 3> 1.2. 2 Derivation of C* Guidance Equation 


The velocity-to-be-gained is defined as the difference between the 
required velocity and the current velocity, i.e.. 


V f *V g (r ,t)-V 

Thus, the derivative of equation 2.3.19 is 


(2.3.19) 


( 2 . 3 * 20 ) 

Now, since Vg, is a function of both position and time, its derivative is 


... , ,n d\/ a dr d\ f q 

V„ (r,i)= — - — + ~ 

~ 3C dt U 


(2.3.21) 


so that equation 2.3-20 becomes 


d'lo dv 

L /„ = V t — 

"? dr ~ dt 


d h-v 


( 2 . 3 - 22 ) 


where V has been used for . 

dt 

The total acceleration of the vehicle is given by V. This function 
is the sum of gravitational and thrust acceleration, i.e., 


Substituting equations 2.3-19 and 2.3*23 an equation 2.3.22 yields 


(2.3.23) 




(2.3.24) 


v, 

V a ~ Yo ~ 


2r ■? 


(2.3.25) 


This equation can be simplified if the definition of required velocity 
is employed. Consider two vehicles, A and B, at the same point in space. 

Both vehicles are to satisfy the same terminal conditions at the same time; 
however, vehicle A has already acquired the necessary velocity in order to 
terminate correctly and is presently in free flight, having terminated thrust. 
Vehicle B, on the other hand, has not acquired its required velocity and is 
still thrusting. Since vehicle A is in free flight and experiences no 
thrust, the following is true: 
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r 



and equation 2.3*25 for vehicle A becomes 


or 


O = 



dt 


"£a 


(2.3.2b) 


2 * ar A -** 


dy 

9t 




( 2 . 3 * 27 ) 


This result is expected since it states that the only acceleration that 
vehicle A experiences is that of gravity. However, since vehicle A is at 
the same position as vehicle B, and since both have the same terminal 
constraints, they must have the same required velocity vectors. Furthermore, 
both vehicles are experiencing the same gravitational field so that each 
point on the thrust trajectory can be compared to a free flight vehicle that 
has exactly the same required velocity and gravity vectors. Of course, each 
point must be compared to a different free flight trajectory since the 
required velocity is continuously changing. 
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Finally, if all points on the powered trajectory instantaneously satisfy the 
required velocity and gravity vectors of some free flight vehicle at the 
same point, then equation 2.3-27 is also true for the powered flight, i.e.. 


w av 

4* J 7 * tt (2 - j - 28) 

along both a powered and free flight trajectory. In this light, equation 
2.3-25 becomes 


or 


dY* 

V ~ V -a 

*1 dr 'f ~ T 


(2.3.29) 


-? t l7-? = ' 4 T (2.3.30) 

d)i„ 

The matrix -jf evaluated along the nominal trajectory is called the C* or Q 
matrix. Thus, in conventional notation 



(2.3.31) 


This equation provides a scheme for computing the velocity-to-be-gained 
from the measurement of the thrust acceleration. The object of the steering 
policy will now be to use the velocity-to-be-gained and its time rate at 
various epochs to drive Vg to zero. The choice of the steering policy 
depends on the quantity that is desired to be optimized} thus, families 
of logics can be proposed. Steering will be discussed in section 2. 3. 1.3. 

In fact, an optimal program including a discussion of the closed loop will 
be presented in section 2.3.1. 3*2.1. 


2.3 .1.3 Steering 

Guidance schemes, in general, and those that use the required velocity 
concept in particular are usually divided into two phases, the atmospheric 
phase and the vacuum phase. The reason for this division is that the 
vehicle cannot tolerate excessive structural loads due to aerodynamic effects. 
Thus, since the major strengths of the vehicle are axial, the loads must be 
near axial during atmospheric flight (i.e., if there is steering, it must 
be "gentle" in nature.) However, after the vehicle leaves the atmosphere 
(or more correctly, after the aerodynamic loads are reduced below a specific 
level), it may be subjected to the more violent maneuvering that may be 
commanded by the guidance system. For this reason, the steering is usually 
run in an open loop manner during the atmospheric phase in order to prevent 
any violent commands that might occur during this phase. The open loop 
steering is, in general, designed to keep the vehicle as close to the near 
optimum nominal path as possible. The nature of perturbation guidance 
schemes requires this closeness to the nominal in order to assure the 
accuracy during the steering phase to be guaranteed. 
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Acknowledgement is given to D. F. McAllister, D. R. Grier, and 
J. T. Wagner whose work on the optimum steering for the powered phases of 
the Apollo Mission (reference 3-5) was used in the preparation of this 
section. 


2. 3. 1-3.1 Atmospheric Phase 

The atmospheric portion of the flight usually consists of a vertical 
rise for a prescribed period followed by a transition turn. The object of 
the transition turn is to rotate the vehicles attitude and velocity vector 
by an amount known as the "kick 11 angle. This kick angle defines a zero 
lift (gravity turn) trajectory through the atmosphere such that near nominal 
(optimal) conditions are attained upon entry into the vacuum phase. One of 
the simplest ways of implementing the atmospheric phase of steering is to 
mechanize a program of vehicle attitude or attitude rate based on sensed 
data. During the kick, yaw attitude is usually held to zero, while pitch 
attitude is commanded. 

A more complicated loop steering scheme that is used during the atmos- 
pheric phase is called velocity steering. In this scheme, the desired 
vertical velocity is 'written in terms of position and time for the atmospheric 
phase of the flight. The true vertical velocity is then compared to the 
desired value and a pitch perturbation command that is proportional to the 
difference is generated. The pitch perturbation command is then added to 
the preflight nominal pitch command (which is also a function of position 
and time) to generate the total command. The primary advantage to the 
velocity steering method is that trajectory perturbations, as reflected in 
the variations with respect to the nominal trajectory, are greatly reduced 
during the atmospheric phase. On the other hand, this method introduces 
smaller stability margins at vehicle vibration frequencies because of the 
more active attitude control system. 

2.3 .1.3 -2 Vacuum Phase 

Once the vehicle is out of the atmosphere, structural constraints can 
be relaxed and the steering system can begin to perform its primary function, 
that of reducing the velocity-to-be-gained vector to zero. One obvious 
method of driving Vg to zero is by thrusting in the direction of V . This 
method could be used for either Delta Guidance or C* Guidance since both 
methods calculate Vg either directly or indirectly. A simple mechanization 
of this steering law is obtained if the vehicle is given an attitude rate 
command that is proportional to a^ x V , where ay is the thrust acceleration 
vector. The verification of aj x 1L sneering can be made by considering the 
rate commands that would be generated for various orientations of ay and Vg. 
The vector ay x V_ is zero if ay and V are aligned and no attitude change 
is commanded. If^they are not aligned^ however, an attitude rate command 
in the ay X V direction is given > and the vectors begin to realign. Since 
the C* Guidan§e Scheme generates V as well as V , another type of steering 
that nulls the cross product of these two vectors is # suggested. In other 
words, Vg can also be driven to zero by keeping the Vg vector antiparallel 
to Vg instead of ay. As a matter of fact, it is shown in section 2. 3.1. 3-2.1 
that the V, x ^ method of steering is optimum for C* Guidance in the sense 
that a minimum fenount of fuel is consumed. 
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A comparison of the steering schemes indicates that these two methods 
in general are not equivalent. This can easily be seen by considering the 
expressions for aj. and V g in each case. For the x Vg case 

£t 

-r y f 


or 




h 

v. 


Employing the C* Guidance equation (2.3.31) it is seen that 
Similarly, for the V x Vg case 



So in general, the steering schemes are not equivalent. The following 
sketches give a pictorial interpretation of the two steering schemes . 
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In each of these sketches the nulled condition is assumed. The appropriate 
steering can be verified if misalignments are considered. 

An interesting adaptation of V_ x V steering to Delt a ^ Guidance by a 
numerical differentiation of the velo city-to-be-gained for Vg could be 
conjectured. (it is recalled that a.p x Vg Steering could be applied 
readily since was not needed). Such a differentiation process would 
introduce a noisy V_, however, and may deteriorate the end performance in 
spite of the more optimum steering law. Simulation studies would indicate 
the best method to be used for a particular mission. 

Other steering schemes using the velocity-to-be-gained require a 
knowledge of nominal steering commands for the entire flight. It should 
be noted that cross product steering made no such requirements. Vg can be 
driven to zero by interpreting its components as error signals. A suitable 
attitude command would consist of the sum of the nominal command and some 
combination of the appropriate component of (typically proportional plus 
integral). The following sketch shows how such a scheme could be used to 
null one component of V„ for either C* or Delta Guidance. 



If a vehicle wanders from the nominal trajectory by a substantial amount, 
the steering commands could be larger than tolerable. This results in 
nonoptimum use of propellant. In order to alleviate this problem, a 
nominal velocity-to-be-gained can be precalculated and a less stringent 
error signal such as the deviation from nominal of V g can be used. Further, 
by introducing a weighting technique, the error signal can be weighted 
lighter at the beginning of the flight when (required velocity) is not 
known so accurately (as in Delta Guidance). 
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2. 3*1. 3*2.1 Optimum Steering for C~* Guidance . In section 2. 3*1*2. 2 the 
C guidance equation was shown to be 

-5 * C * ^5 = ~£t (2. 3*3*3) 

The mechanization of this equation with a steering law constitutes a closed 
loop guidance scheme as illustrated in the following sketch: 


c * - c 

* Ct) 

J 




The object of the following analysis will be to determine a steering policy 
for this loop. The simplest steering policy would be to drive the thrust 
acceleration in the direction of the velocity-to-be-gained . Usually, however, 
it is desired to maximize the performance in some sense. The following 
derivation presents the formulation of conditions necessary for a steering 
scheme to be optimum and shows that cross product steering ensures the 
reduction of the velocits^-to-be-gained vector to zero while maximizing the 
burnout mass of the spacecraft (minimize propellant consumed). It should be 
noted that such an optimization problem requires the use of Pontryagin's 
Maximum Principle. Since it is beyond the scope of this monograph to present 
this principle in rigorous detail, the reader is referred to reference 3*23 
for an introductory explanation. 

There are many ways of formulating the state equations for the variable 
mass vehicle. However, the ideal approach involves the use of variables 
which uncouple the equations being processed. Appendix A presents several 
approaches designed to accomplish this objective and substantiates the 
choice of variables used for the following analysis. Let 


*/ -- 



(2.3.34a) 

Xa - 



(2.3.34b) 

= 



(2.3.34c) 

X* = 

<3 I iP | 

= V e JUl f-P^) 1 

(2*3*34d) 
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(2.3.34e) 


U, 


LL^ 


l Tz. 


*T 


ex- 


it 


* - 


where 


_ ImCtQl 

i*»CO 


^ U) 


( 2 . 3 . 34f) 
(2.3. 34g) 
3 ( m < 0 ) -34h) 


rnCt) 

= components of 
G-rj, = components of a 7 
^ t = |&t| ~ thrust magnitude 

Ve. = escape velocity of exhaust gas 

It is desirable to express equation 2.3-34 in principle coordinates, i.e., 
eigenvector directions of C*-, in order to simplify the algebra involved in 
the analysis. If such a coordinate system is employed, equation 2*3-3 -3 
becomes 
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A, 

o 

o 


o 
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o 
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■ 



a T, 



— 

(Xr z 


_ v 53_ 


0-t 3 

— — 


or 


X, O 

O A 7 


0 

o 


(2.3-35) 


(2-3-36) 


where all components are taken to be along principle axes. 
Now, Note that 




m(o) mCi) 
too (o) 


~ Ve. 


rttt) 

*^Ct) 


(2-3-37) 
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But, since * v '(t)is the thrust of the rocket, i.e., 

u, ”C±2 = -ar 


(2.3.3s) 




and since 


met) 


= u. 


X* = Ve 

Thus, equation 2 . 3*34 can be written in the following state variable form. 

x, - h,X, + c/* 0f t (2.3.40a, 

X z - A 2 X 2 * i> e ,U^U 2 (2.3.40b 

X 3 s A 3 x 3 y- u 3 (2.3.40c 

X* - (2.3 -40d 


or in matrix notation as 

X ( A, o o o 

X 2 o A 2 o o 

X 2 o o A 3 0 

X* 000O 


X, 000 

X 2 o o o 

X 3 <2 t/e^ o 


X = AX +8 1/ 

Finally, the adjoint equations to this set ate 

^3 A/ o 00 

A _ _ 0 A 0 0 

A o o A 3 o 


0000 


£ * - A r /? 


See section 2 . 3. 2 . 1.3 for a discussion of adjoint equations. 
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At this point the generalized Hamiltonian defined below is introduced 
H = P T - X * P, X, + P 2 + P 3 K 3 + P+ X„ 

~~p,iKt, + v & u ¥ u,) + r 2 (A 2 X Z + i/ & ^u 2 ) (2,3,45) 

^CA 3 X 3 + Ve.L/+U 3 ) + 

and the performance function is defined to be the burnout mass, i.e.. 


S T = X+CT) = 14 jUl 


rn Crj| 
[/>7COjJ 


(2.3.46) 


In standard functional notation. Equation 2.3.46 is 


S T * C -X T =[00 0/ J [*' T 


Z T 

<tf "T 


(2.3.47) 


Now, a necessary and sufficient condition for maximizing the functional, S,p, 
is that the Hamiltonian, H, be a minimum with respect to the control vector 
at every point of the path subject to the terminal condition that 


pj * £(t‘T) = - C T 


(2.3 .48) 


Thus 


T) = 'C 7 * = 


O 

o 

o 

-/ 


or 


P+Ct = T) - “/ 


(2.3.49) 


(2.3.50) 


But, since P^ = 0, then P^ = -1 for all t so that equation 2.3-45 can be 
rewritten as 


H = A, P, X, + P z X z t /\ 3 P 3 X 3 t V t U r </, v /J U z *P 3 U 3 -/] 

4 P - * 


4,o 0 


X, 

(r 


0 0 


X* 

+ P 3 J 

“z 

0 0 Aj 


*3 

L 



(2.3.51) 


- i 
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or 


H = P* T C~C*) X* f l/e U+\E* T (d* -/] (2.3.52) 


where 



Cl t 

/^rj 


This is the H function that is to be minimized with respect to the control vec- 
tor at every point of the path. Before proceeding, it must be remembered that 
u^, which is defined to be m /m, is limited to the allowable range of 


-D < (2.3.53) 

where F 1 

° ] 

This constraint states that there is a finite amount of thrust available and 
that the direction of thrust cannot be reversed. In order to determine the 
control that minimizes H, it is necessary to investigate alj. possible values 
of the second term in^he above equation : i.e. , V e u ^ P 5 ' - 1 . Consider 

the two cases, (l) |P*'| >1 and (2) 0 4|P*| 41: 

(1) lE l > 1 • For this case, two choices of control are possible. One 
choice of control would be ^ ~ 0. Under this choice the control term of 
equation (2.3*52) would be zero. A smaller value of H can be realized, how- 
ever, if is chosen to be -D and [P* Tu* _ qj ± s chosen to be some positive 
quantity. Under this choice of control 

(2.3.54) 


w. m ->t 

is most negative if p" i u - 1 is made as large as possible; i.e., in the 
limit 

-l£'l or u* - ^ 

Hence, if P*= 1 the control that minimizes H is 

P * 

~ U* = JjhTj (2.3.55) 

(2) 0 < |P'/ < 1 . For this cgse, since 0 c /P"/<1, then fP* - lj 

must be a negative quantity (P' r ^ u" < l). Thus, cannot be positive and the 
minimum value of H is attained for the choice = 0, i.e., if 0 < JP"| <1, the 
control that minimizes H is = 0. 

In summary, the optimal control law is to burn at maximum thrust in the 
£* direction as long as Jpf/ >1. When /P'7 < 1, terminate thrust. Thus, the 
thrust vector during powered flight may be written as 
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(2.3.56) 


P * 



where 

= l*rl 

It is desired to determine the optimum steering in terms of the velocity- 
to-be-gained. Thus far, it has been determined in terms of the adjoint varia- 
bles. In order to accomplish this transformation, it is first necessary to 
write the solution to the adjoint equations. 


p, 

- o< * 

(2.3.57a) 

Pz 

* /3 <2f A z^ 

(2.3.57b) 

P 3 

= 

(2.3.57c) 

z ^ p * | = j^ 2 e‘ 2A,t t/3 2 & 2Azt + y 2 c" 2 * 3t ] ^ 

(2.3. 57d) 

Now, the state equations (2.3.35) 

can be written as 


r 

a / v 3/ “ a T Y - 

(2.3.58a) 

v %z * 

Az Vcj 2 ~ <X T 

(2.3.58b) 

= 

?3 

A 3 v^ 3 - a T - 

(2.3.58c) 


Using results of Section 2. 3. 2. 1.3, these equations can be integrated from 
t = 0 to t to yield 





Ct-nr) 


a-r(^) 

nr) 


P/C'T) ( 2 *3. 59a) 


V«J 2 C a - ^crj (2 * 3 * 59b) 
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v<3 3 (0- v^o) ^3 i -J r e .^~ r ) P 3 cr)dr (2.3.59c) 

o 

This set of equations can be evaluated from time t = 0 to the cutoff time, T, 
by substituting the terminal value of the velocity— to— be— gained (zero), i.e., 

V T) = 0 V 01 ^ ~ 3z£21r' (r)</r (2.3.60a) 


or 


l/^zCT) -0 - l/^Co) e.^ T - f e AzC7 ' T) &z^21 p cr)c/y ( 2 . 3 . 60 b) 

t 


^ 3 (r) - vy Co) e" 1 '-/ 


r^AJ-7-nr) OrC'y) 


±Cp) 


P x CP)dy (2.3.60c) 


e A ' T i * y /e A,c7 "' r; e A '' r a"r ( 2 . 3 . 62 a) 

O^Co) e.^-iSj' e ' # *-C T ~' r - ) Q ^ C ^j a (2.3.62b) 


v^Co) e A > T = yy'e 


ZUT-T) a T cr) „-4,r 


zC'r) 


^ (2.3.62c) 


Now the constants at t /3 t and 2" can be found as soon as the integral 
in equations ( 2 . 3 . 62 ) is evaluated; and therefore, equation 2.3.56 can be 
used to find the thrust policy as a function of the velocity-to-be-gained. In 
general, the integrals of equation (2.3.62) are not easily evaluated. However, 
the exponential can be expressed as a Taylor series 


- Z\,^_ 


= / - zA,^r + 


*A,T‘ 


Z{ 


(2.3.63) 
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and. higher-order terms can be neglected if 'X 4 4/ . This approximation 
is not very unrealistic -when the burn time is compared to the time constants 
of the adjoint equations. Under this assumption, the first two terms of the 
Taylor series suffices in equation (2,3.62), 

Ckj C/- ZA,r) dT = Co. ) (2.3.64a) 

* ~r cc c r) 

O - * h t) -c/r = l/^Coj (2.3.64b) 




Q t 

-aCr) 


dr = 




(2.3.64c) 


Since the method of integration of equations (2.3.64) is the same, only 
one integration will be performed. The solution to the others mil be written 
by analogy. This integration mil be performed by parts after first writing 
the equation in the' following form: 


o( 



T OrCr) 

ZCT) 


cir 


-fzA.T 


a r cr) 

icr) 


dr 


^ l /$Co) (2.3.65) 




T a~<~r ) 

kir) 


dr - zA, 


< ^61}dr 

*cr) 


-V^Co) (2.3.66) 


and substituting 
to obtain 

(A \ XCT) — Z A, 


ICT) - f 


Wo 
Or £T) 

H cr) 


dr 


t tCt) -Jjcr) dr 1 

t? J 


r x d^,Co) (2.3.67) 


Now, assuming that l(T) can be approximated linearly over the interval 0 < t 4 T 
as 


J I CT)d'r = (2.3.68) 
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(This again is not an unrealistic approximation since the burn time has been 
assumed small compared to l/A.) allows equation (2.3.67) to be written as 


CK 


{ 


I(T) -2A.T-ICT) + 


i 




Co) 


or 


C,[l(T) — A,T’I(T) j ^ * /- A,tJ Z(T)= Vf/o) 

A similar technique can be used for /3 , and ?f. The results 
below. 


<*[/ 

-V 

Jia) = 

,(0^ 

4- 

- a 2 t 

]xcr) = Vj 2 

.Co) 


a 3 t] 

II 

K' 

V — / 

H 

(°) 

<X , £ . 

, and ' 

7 are 


dv 




= (T- 

-\T)Z(T) 


/2 - 




& 

= 77 

- x z T)rtT) 


X - 




0 

<7- 

A 3 T)X(T) 



But, since the approximations xT«l have been made 


/ 


/- AT 

4 


= !+\J 


is valid, and equations (2.3.72) become 


(2.3.69) 

(2.3.70) 
are stated 

(2.3.71a) 

(2.3.71b) 

(2.3.71c) 

(2.3.72a) 

(2.3.72b) 

(2.3.72c) 
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G< 

0 

V 


-C/*A,t) 

c ° J 

(2.3.73a) 


Z LT) 



v^c°) 

(2.3.73b) 

= (// A*r) 

I Ct) 

- Lit A 2 t) - 

i cry 

(2.3.73c) 


Now that the coefficients to he adjoint equations have been determined in 
terms of the initial velocity -'to-be-gained and the eigenvalues of the C* 
matrix, the initial unit control vector u*(o) can be determined by the defini- 
tion in equation (2.3.55b) as a unit vector in the Pf direction. 


cm = 


/ 

ICT) 2 CO) 


V^co) 
(/*A z t) V$z C o) 
VtA 3 T) V$ 3 (o) 


(2.3.74) 


f 

(It \ t) 

o 

0 

- 

ty \o)~ unit 

0 

(!+ Azj) 

0 

y^(o) - 

1 

o 

0 

(i+A 3 t) 

4 


However, since u* is a unit vector, it is not necessary to include the coeffi- 
cient izi5) • Rather, it is sufficient to write u*’ as 


(2.3.75) 


with the understanding that once the vector 

(t*A,T) v^Co) 

V$ z (0) 

(/*\ 3 r) V$ 3 co) 

is determined, it will be scaled to a unit vector by dividing by its magnitude. 
For subsequent times, the optimum control is defined by the condition 


U*Li) = unit , 






(KM 


which, under the substitution of the time-to-go (t ) defined by t gQ = T-t 
becomes 


6/ ft) - Uni t ■) 


go 

(ltA,t $0 ) 0 0 

o C/tAz'C-So) o 

o o (a Ajtjo) 


tfi(t) 


(2.3.76) 


89 


Equation (2.3.76) gives the optimum thrust vector orientation to drive Vg to 
zero at time T while minimizing propel] .ant consumption for a C guidance 
scheme. 


It is now desired to find a steering scheme (no attempt will be made here 
to show that the resultant scheme is unique) that satisfies the condition of 
equation (2.3.76). This objective will be accomplished by assuming a candi- 
date steering law for C* guidance based on nulling the vector product 


x V 

-s- % 


and proving that this steering satisfies the conditions of equation (2.3.76) 
(thus, is optimum in the sense that it minimizes the propellant consumption) . 
Nulling the quantity Vg x Tr in effect forces V to be anti-parallel to V for 
all t, 0 i=t ^tg Q . This behavior is due to the ract that the cross product 
is zero if the two vectors are parallel, anti-parallel, or if one or both 
vectors are zero. While thrusting, the only possibility for nulling this 
cross product is if Vg is anti-parallel to Vg. That is, this steering policy 
produces a thrust vector in the Vg direction which is anti-parallel to a unit 
vector in the Vg direction (for all t) ; i.e. , 


or 


V % 



(2.3.77) 


The desired thrust can now be found by substituting equation (2.3.77) into 
equation (2.3.33) to obtain 



+ C* 



(2.3.78) 


or 


— r -- + 






(2.3.79) 


But, in principal coordinates, the C~ matrix is 


-C 


* _ 


A, 

o 

o 


- 

o o 
A* o 
o A 3 


(2.3.80) 
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Thus, equation (2.3 *79) becomes 



(2.3. SI) 


Finally, 

( 2 . 3 . 81 ) 


since the time-to-go can be approximated as 
can be -written as “ 



I+A^o o O 

o o V 2 

o o I + 


zL„ , equation 


( 2 . 3 . 82 ) 


and the unit vector in the thrust direction is 

u*=unit‘ O I + o 

Q O J + 

J- 

But, this result is exactly the equation that must be satisfied for ''optimal” 
performance as specified by equation (2.3.76). Thus, the cross-product steer- 
ing scheme is optimal in the sense of minimum fuel consumption for C" guidance. 

The use of ILxV to generate a proportional steering rate command suf- 
fers from two related §ources of error. First, as V g approaches zero, the 
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cross product also approaches zero, thus reducing the gain of the control 
loop. Second, since Vg is not exactly coincident with the acceleration vector, 
a small turning rate must be given to the acceleration. vector in order to 
drive V to be anti-parallel to Vg. This turning rate can be calculated by 
considering the time derivative or a unit vector in the V_ direction, 1 V£T . 

CD 9 "~ V O 



(2.3.84) 


Since l V g is a unit vector, its magnitude does not change in time and its 
total derivative is 



(M X 


(2.3.85) 


Further, if is taken as the pitch and yaw rotation rates of the unit vector, 
then u> must be perpendicular to l V g. In this light 

~ = ~ *[- X ^3] (2.3.86) 


So 


& = 


5 Vi 







LO - 


/ j/j 


v i 


(2.3.87) 


Hence, the turning rate of the vehicle is proportional to the misalignment of 
Vg and Vg which is the policy of cross-product steering. 

Since equation (2.3.87) prescribes infinite turning rates as V -► 0, it 
is customary to command the turning rates as 


±0 -- K ( y 5 X ) 


( 2 . 3 . 88 ) 


where K may or may not be a function of Vg. Writing uj in pitch and yaw com- 
ponents, the command rates become 

OOp - ~K x t>, ft * VjK ) (2.3.89a) 
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(2.3.89b) 


^ = 


K ( u 3x ~ 


In summary, the optimum steering for the velocity-to-be-gained C* guidance 
has been found to be Yg x L (cross product) steering. First, the steering 
conditions for minimum fue3?c on sumption -were determined. Then, Vg x V steer- 
ing was shown to satisfy these requirements. Finally, it was shown thft 
vehicle-command rates must be proportional to Yg x Vg in order to maintain the 
anti-parallelism of V and V_ as required for cross-product steering. 

o o 

2.3.2 Linearized Perturbation Guidance 


Linearized Perturbation Guidance is a scheme that takes advantage of the 
assumed fact that if a reference trajectory and control is specified for a 
flight, the actual trajectory will be "close" to the nominal trajectory. This 
"closeness" enables the nominal-nonlinear equations to be reduced to a linear 
set of first-order differential equations with time-varying coefficients. The 
techniques of linear analysis can then be used in order to relate control 
deviations to position deviations from nominal as the vehicle travels along 
the trajectory. The section begins with the basic linearization technqiues 
and presents the solutions to the linear time-varying second-order differen- 
tial equation. The solutions are then related to the guidance problem. 

2. 3. 2.1 Formulation of Linearized Perturbation Equations 

2. 3. 2. 1.1 Linearization Technique . The nonlinear equations of motion for 
a vehicle can be written as 


X = F(i,X(i), 1/(0) 


(2.3.90) 


where 


X(t) is the state vector 
U(t) is the control vector 
t is time 


F is 


a function vector 



It is assumed that a numerical solution has been generated with some prespeci- 
fied control law that is optimum in some sense, and that this nominal trajec- 
tory and the corresponding control as functions of time are 

%*(t) 

y*ct> 
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During the actual flight, the vehicle will follow a trajectory X(t) which 
differs from X*(t) due to errors in the system and in the model used to gener- 
ate the nominal solution. Thus, the associated control, U(t), required to 
correct for these errors must be computed. The perturbations in the state and 
control vectors are 


6ZU) = X(t) ~X*(t) 

(2.3.91a) 

8y(t) = uct)-u*ct) 

(2.3.92b) 

So 


dt 1 " J dt dt 

(2.3.92a) 

and 



(2.3.92b) 

Now, the nominal trajectory and control relationship is 


- X*Ct) = F(tJC*U), ir*(D) 

dt - 

(2.3.93) 


Thus, from equation (2.3.92a) 


— = F(t,X(t>,ycO)-F(t^*(i), (2.3.94) 

dt 


However, since the actual trajectory is assumed to be close to the nominal tra- 
jectory, the term F(i,X(i), If (t)) can be expanded in a first-order truncated 
Taylor series about the nominal; i.e., 

* * 


6X CtU 


d£ 

HJ 


6UCi) (2.3.95) 


where 

exit) = xct) - x (« 

51JCt)- l/U)-U*U) 
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d£ 

dX 

d£ \ 


the Jacobian matrix of F with respect to X 
evaluated along the nominal trajectory 

the Jacobian matrix of F with respect to U 
evaluated at the nominal control law 


Note that terms for deviation in time are not included. This deletion is the 
result of the fact that a fixed final time is assumed. 

Equation (2.3.94) now becomes 




* 


+ 


av 


i'xrct) -£djX*a)\r*a)) 

y 


(2.3.96) 


or 




bX 


sxa) + 

* 


iv 


sva) 

* 


(2.3.97) 


At this point, a rotational change will be accomplished by substituting 


Ad) 

B(t) 


ax 

af 

Tx 


%(t) 


SX(t) 
S VCt) 


so that equation (2.3.97) assumes the familiar time-varying linear differen- 
tial perturbation equation 


i.(t) = Ad) XU) * <0(0 u(t) (2.3.98) 


2. 3. 2. 1.2 Homogeneous Solution . The homogeneous solution to equation 
(2.3.98) is 




(2.3.99) 
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■where Xq is the state vector at time t = t Q and 0 (t, t 0 ) is the state-trans- 
ition matrix -which is defined as the solution to the following matrix differ- 
ential equation 


J - A{t) J (2.3. ICO) 


Equation (2.3.99) can easily be seen as the homogeneous solution by taking the 
derivative of equation (2.3.99) and employing the definition of 0(t, t 0 ) 


=A(t)§( tjtJXo 

' A Hi K(.ti 


( 2 . 3 . 101 ) 


So 

*60 -- A(t) X(t\ (2.3.102) 


which is the homogeneous equation. 

Two of the properties of the state-transition matrix will be useful later 
in the development. They are 


‘<$(*3 A) • (2.3.103) 

- J identity ma-ir/'x.) (2.3.104) 


2. 3. 2. 1.3 General Solution . The state deviations can be controlled by 

proper choice of the control deviation, u(t) . However, in order to determine 
the effect of the control on the state at a later time, it is necessary to 
know the general solution to equation (2.3.98). As will be apparent, the 
general solution to equation (2.3.98) alone is not sufficient knowledge for 
determining a control law in itself. 

The general solution to equation (2.3.98) is obtained by the introduction 
of another set of equations called the adjoint equations. These equations 
reduce the problem to a straightforward integration problem in which there is 
no cross coupling in the state variables. The proper choice of the adjoint 
equations can be determined by considering the derivative of the product 
_aI ( t ) , where (t, t Q ) is a matrix of variables for the 

adjoint system. (This matrix will be defined during the steps which follow. 

At present, the elements of this matrix are unknown.), and where "T" represents the 
(conjugate) transpose of a matrix . 
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That is, 

^[aUuxw] 5 

= X^i^lUt) ^A(t,i D )A(i)X(i) + A.(t,i,)B(i)U(i)( 2.3.105) 


Now, since the elements of A(t, t Q ) have not been specified, they can be 
selected to satisfy the equation 

K{ tjt 0 )= - X.(ijt 0 ) A(t) j 

( 2 . 3 . 106 ) 


with the boundary condi tions_^y_(t 0 , t Q ) = 1 
Thus, equation (2.3.105) reduces to 


y.t) 


AitJo) 3(t)U(t) 


(2.3.107) 


The integration of this equation is now seen to be straightforward. The 
result is 


* jA(% > to)8(%)M(S J ^1(2.3.108) 


Now, employing the fact that A (t Q , t Q ) = I, the solution for x(t) is 
obtained as 





(2.3.109) 


Hence, the general solution for x(t) for all t > t 0 has been determined in 
terms of the adjoint matrix. 

Now, consider the derivative of the product 

AdjJ $(t,0 . 
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That is, 

d_ 

dt 


{X ( 4,0 §(*,*.)]= < AttjU) §UAJ (2.3.110) 


But, substituting equations (2.3.100) and (2.3.106) reduces this equation to 


i,iJ$(4,o] = -XUjtJAU) $(*,0 (2.3.1U) 

+ KLtj t Q ) A(t)<§>{tjt 0 ) -O (A/ulf natr/x) 


so that 


*XU 0 >+o) $ (Ujto)+J Odt -- 1-1=1 (2.3.112) 


But, if 


then 


or 


A r (^0 $(*,£„) -I 


$(t,4.) - -K ( tjQ 


t a ) = ^7w.) 


(2.3.113) 


(2.3.114) 


(2.3.115) 


Now using the relationship between 0(t, t Q ) and _A.(t, t Q ), equation (2.3.109) 
can be written in terms of 0(t, t Q ) as 

yjt) = (2.3.116) 

■t-o 
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(2.3.117) 



But 


AUi.) = <£"(^0 (2.3.118) 


So, finally 


X(i) = SlUjic 

't 

/ $ (i,S)0(S)U($)c/$ 

JL 

4 


(2.3.119) 


Equation (2.3.119) is the general solution to the differential equation 

X 6t) = A U) a C-i) (2.3.120) 


Note that this solution is independent of the adjoint parameters and requires 
no inversions. 

2. 3. 2. 1.4 Fundamental Guidance Equation . The derivation in Section 

2. 3. 2. 1.3 gave the general solution to the first-order perturbation equation. 

A slightly modified form of this equation, known as the fundamental guidance 
equation, is used extensively for performing error analyses, generating require- 
ments, commanding corrective maneuvers, making linear prediction, etc. This 
section will derive this equation and give some of the interpretations of it. 

The general first-order perturbation equations will be rewritten for con- 
venience. 

X ftJ =/!C*) X C-t) 1-3U)CS U) (2.3.121) 

eft 

Now, since 

A (ijto) - - a t (i) a(v.) 

an adjoint system can be defined as 

X(t) = ~A T U) Xct) (2.3.122) 


where X (t) is the state vector for the adjoint system. 

Thus, an analysis similar to that used in Section 2. 3. 2. 1.3 can be used to 
derive the useful combination of the two systems of equations. 

Consider the derivative of the scalar product of X. and x. 
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=<*,Z>+<d,2£> 

= <-/4 T A ,x)t(A ; /!X+^) 

- <^~A r X> -/-(A >*)(} + (* >BU ) 

= < *: J&y > ~ <* r 6;“} (2.3.123) 

where ( A > % > = X ' x r inner product, of two vectors 

Integration of equation (2.3.123) from some initial time, t = t 0 , to a final 
time, t = T, yields 

X r (T) x. CT) - X T (ij £ U ft ) » ^ X 7 (r> B(r)y(T)dr 

or 

X t ct)^ct>= X T ( i j ;aO + / Vf r) BtT)uc*)d* (2.3.124) 


This equation is generally denoted by the fundamental guidance equation. 


So far, no constraints on the systems of equations have been specified. 
However, the proper specification of these constraints reduces equation 
(2.3.124) to a more useful form. One constraint that can be specified is the 
terminal state deviation, x(T). This constraint can be expressed as some max- 
imum allowable deviation at the end of the powered trajectory. It is noted 
that though x(T) does not explicitly appear in equation (2.3.124), it does 
appear in the form of the inner product of <XCD, x(.t)> . Thus, if the appro- 
priate choices of the components of A (t) are made, a particular component of 
x(t) can result from this inner product. For instance, if the inner product 
is desired to be X]_(t), the choice of A (T) is 


A(7) = 


O 

O 


L OJ 

It should be apparent that in order for x.he left side of equation (2.3.124) to 
be written explicitly in terms of x(T), a different adjoint vector is needed 
for each component of the state vector. 
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rz, ct> 

xl(n^cr)= [/ 00 ... 0 ] x z ct) 


x„ cn 


= 2,(73 


[X<7r 

(T) x<r; = [ 0 / 00 ... oj ^< 7 ; 


cn 


-X^7)J 


r z, f7>- 

£ n (T)x(T) = [ 000 ... 1 ] ^(r) = 2 „cr; 

.X n (T)_ 


Equation (2.3.124) can be written for each adjoint vector that was selected in 
these constraint equations. Written in matrix form 

X(T) X (T) = Xuj z(t a ) +fX(T) BCT)U(T) dr (2.3.125) 
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A r (7> 


AiO~- 


Ait) = 


X W 


V T) 

= 


L 

Auj 

r 

VO 

= 

An'iO 


X a) 

r 

i z u) 

= 

A r it). 

— n 

L 


too... 01 


= 1 


ooo . . . / J 


VO X (2 U 0 )...A ln (0' 
X 2< (£„) A.^^) ■ • • X 2n C t e ) 


'•‘■O' nz^ c ® nn * 

X„«) X 1Jt (0.. ■ x,„(t) 

x 2 ,(i) x„rt)...x 2 „«) 


Now that the boundary conditions of the adjoint equations have been spec- 
ified, these equations can be integrated in reverse time from t = T. 


A(t)= -A T (t)AU) 


(2.3.126) 


A < T ) = I 


(2.3.127) 


Thus, for a given state deviation at an arbitrary time, t, the only unknown 
term in equation (2.3.125) is the integral. (All other terms are known since 
the terminal and present state deviations are specified and the adjoint equa- 
tions have been integrated.) Hence, the control that is necessary for a pre- 
scribed terminal- state deviation is implicitly contained in the integral of 
equation (2.3.125). This information is useful in determining possible control 
for the remaining portion of the powered flight. 

A simple example would be the determination of a constant control vector 
(for the remaining powered flight time) that would satisfy the terminal state 
deviation specification. Consider the following: 

f _A(T)B(T)U dr = D 

where D is known from the specification of present and terminal-state 
deviations and the adjoint equations. 
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Since u was arbitrarily selected to be a constant vector for this simple example, 
it can be factored out of the integral 


The integral can now be 
control vector is thus 


[JmtI 8(T) c/rj U = D 

found since A (t) and B(T) are known. 


*-[A r < 


j\(T)Q(.r) jr 


-/ 

D 


2. 3. 2. 2 Optimum Control 


The constant 


The basic formulation and solution to the linearized perturbation equation 
has been presented in Section 2. 3. 2.1. Although this information is necessary 
for a complete understanding of the problem, it does not explicitly give a 
solution to the control problem. In other words, the material presented so 
far is suitable for use with a control that is known beforehand but is not 
sufficient to determine a control law as such. The theory must be extended 
further in order for it to be useful in control law determination. This exten- 
sion will be realized by optimizing the performance in some sense over the 
entire trajectory. 


The optimum, control problem has been studied extensively in recent years, 
and as a result, a large amount of literature is available on the subject. 
However, since the methods of analysis, while different, are analogous and since 
it is impossible to discuss all of the available material, several of the more 
important variations will be presented. 


Acknowledgement is given to W. F. Denham and A. E. Bryson whose work on 
terminal control for a minimum-mean square deviation from a nominal path 
(Reference 3.19) was used in the preparation of this section. The material 
presented an optimal linear control generally seen in text books on optimal 
control of discrete-time systems. Reference 3.26 presents more advanced con- 
cepts in optimum control theory as well as the basic deviations. 


2. 3. 2. 2.1 Terminal Control for a Minimum-Mean Square Deviation from Nominal . 
In Section 2. 3. 2.1 the linearized perturbation equations for a nonlinear system 
were derived and the general solutions were written in terms of the state trans- 
ition matrix and an integral that involved the control deviation. However, the 
determination of the control deviation that is necessary to correct a state 
deviation is not easily accomplished since the control must be known beforehand 
in order to evaluate the effect of the control on the state deviation at a 
later time. Further, it is usually desired that a control policy be optimum 
(in some sense) and that the on-board computation minimized. Thus, this dis- 
cussion will be initiated with the presentation of a convenient control scheme 
for the linearized perturbation equations. 


In this scheme, it is only necessary to multiply the state deviation by a pre- 
computed matrix in order to find the desired optimum control deviation required 
to guarantee the correct terminal state. However, in order to achieve this 
simplicity, the burden of the computation has been placed on the pre-flight 
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VARIABLES 


simulation. Thus, in addition to the optimum nominal trajectory computation 
that must be performed, it is also necessary to determine the matrix multi- 
plier, -A_(t), as a function of flight time. This latter objective is not 
easily realized since the solution depends on many parameters such as terminal- 
condition constraints, optimization desired, weighting, etc. However, the 
following discussion is designed to present the derivation of the _A-(t) matrix 
for the control such that the mean value of a positive definite quadratic 
form in the control variable deviations is minimum. 

As in any terminal control problem, it is necessary to define the terminal 
conditions that must be satisfied at the end of the control time. If the nomin- 
al trajectory and nominal control are followed for the entire powered flight, 
the Vehicle will arrive at the terminal point with the desired terminal condi- 
tions 

where ^ is a column matrix of "P" known functions of X(t) and t, (P ±/u) 

and T — terminal time. 

Since nominal conditions are not followed exactly, it is necessary to define 
the terminal point by choosing some scalar function fi-[X(T) y T] } which 
satisfies the desired terminal constraint, i.e., control is terminated when Si. 
is satisfied. Such a formulation enables the terminal point to be defined in 
terms of a desired constraint equation rather than in terms of the independent 
variable, t. This feature is advantageous since the nominal-terminal time, T, 
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may not be the time at which CL is satisfied once deviations in state and con- 
trol are experienced at any point in the trajectory. 

Using previous results, the linearized perturbation equations are 

- £<£) = 460 (2.3.128) 


In Section 2. 3. 2. 1.4, the general solution to this equation was found (with 
the aid of the adjoint equations) to be 

XiT) x(.r) -SiUJ Air) B(r) (j(r) dr (2.3.129) 


is a matrix of multiplier functions. 

conditions must now be given to the matrix of functions, yv^(t) . 
accomplished by introducting two different matrices, jC r (t) and 
each of which satisfies the terminal conditions 


where -Aj'(t) 

The boundary 
This step is 
A 'silt) 


A (T) 


A* (T) 



(2.3.130) 


If these conditions are used in equation (2.3.129), the following is obtained 


(6P) 

t*T 



( ) dt t (A" x) 

* *-*. 


(2.3.131) 


where 


(6 Cl) 


t=r 


= di +(A t „x) 

A i 


(2.3.132) 


<S?{. r -<K*Lr 
( * n K- r =< A «*-lr 
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and where 


A 7 '(t) denote the matrix of multipliers defined by 

matching the boundary conditions for both 
V and -Cl , respectively. 

The terminal point is defined to be that point when A is equal to some 
specified value, i.e., the control is terminated when dfl = 0. However, small 
deviations from the nominal trajectory produce changes in the value of T, dT 
and in the total derivatives of V' and A . 

df = (5^)1 + Yc/r (2.3.133) 

1 1= T - 


where 


dSl~~ ( 6X1) 


+ SLdT 


t-T 


r 


ay dx 
3± * ax dt 


A 


3X1 

at 


ML d£ 
3X d t 


t= T 


(2.3.134) 


(2.3.135) 

(2.3.136) 


Substituting equation (2.3.131) into equation (2.3.132) yields 

dr = A^ dt+ fdT *■ ( £ 0 ) _X_(t 0 ) - O (2.3.137) 

r T 

dSL- SXdT + -O (2.3.138) 


It is now desired to have these terminal constraint equations independent of dT. 
This objective is achieved by solving equation (2.3.138) for dT and substituting 
the result into equation (2.3.137). 


r 

df -- j'yZ^e>(i)UCt)dt (6,)x(0 =0 

to, 


(2.3.139) 
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■where 


A yjl. ^ V _Q_ 


is again a matrix of Lagrange multipliers. Equation (2.3.140) determines 
the values of the elements of A 7 ^ , the boundary conditions being given by 


v ,-sJU -X l*k\ 
A ^ ,} Vax A 3 X ) t _ r 


( 2 . 3 . 141 ) 


Now that the constraint equations have been determined, a performance 
index must be specified in order to establish some measure of optimization. 

In general, though it is noted that the performance index need not be limited 
or to the control deviation over the nominal trajectory, these criteria serve 
as very effective measures of the performance of a system. (Any positive def- 
inite quadratic form can be considered to be an allowable term in the perfor- 
mance index as long as it is associated in some way with the performance.) 

For purposes of illustration, a simple quadratic form of the control 
deviation will be used for the performance index. This form is the type which 
would be employed if fuel consumption was of major concern. Assume 


V - j' uct) r Y ui±) dt 


(2.3.142) 


where 7 fs an arbitrary symmetric, non-negative weighting matrix chosen by 
the control engineer and rewrite equation ( 2 . 3 . 139 ) as 

dH' - / \ r rJX _8U)Ua)clt - A^(OX_(0 =° ( 2 . 3 . 143 ) 


Now, variations in V with respect to u can be found once the equations are 
joined by another matrix of multipliers, 

v = f(u T TSu - V r yC yrL 8 u) dt * v T [c/y'-A^OJLCt.)] (2.3.144) 


&V j (BUj X ~ V Si T rs ^B) 6 U^)dt (2.3.145) 
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But, for an extremum in , = 0. This equality occurs -when the integrand of 

equation (2.3.145) vanished, i.e., 


ZU T *'-V T XV_ a £ = o 


( 2 . 3 . 146 ) 


Now equation (2.3.148) can be written 

d ¥ ~f 8 % 8 dt ~ X (6 0 ) — o (2.3.147) 

to 

or 

dY XCO =■ jJ'(to) ^ (2.3.148) 


where ^ 

=/- ' J Zr si e*~ , BA if ,j L clL (2.3.149) 

to 

Now the multiplier matrix which was introduced can be seen to be 

V(±o) = Z J~'Cto)^dr -Ay^to) X( 0 ] (2.3.150) 


Finally, the control deviation can be found in terms of the state deviation 
from the transpose of equation ( 2 . 3 . 146 ). 

dY ~ A’^^XCto)! 5T Cto)l 

J J (2.3.151) 

u = X~'B r a)tyg J- , U.'i\dV 

Equation (2.3.151) leads to the final relation between the control and 
state deviation when d<^ is set equal to zero. 



Ult) = J((U 
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In the way of explanation, d* specifies the terminal constraint deviations 
that are desired based on currently available information. For simplicity, 
these desired deviations ■will not be considered here. The interested reader 
•will find information regarding terminal constraint modifications in Reference 
3.19. 


Equation (2.3.152) can be mechanized in two methods. The first method 
can be thought of as a discrete-time system, in which a major computation cycle 
occurs at intervals that are spaced by minor computation cycles. The second 
method utilizes the computer's speed to run the problem as though it were a 
continuous control problem. The following paragraphs will describe the calcu- 
lations that are necessary to mechanize this guidance scheme in both the con- 
tinuous and discrete modes. 

First, the choice of the no min al trajectory that satisfies terminal con- 
ditions must be made. Then, the elements of the A(t) and B(t) matrices must 
be determined along the entire nominal trajectory. Having determined A(t), 
the matrix of functions ^ (t) can be determined by integrating the adjoint 
equations backward in time from the terminal point with the boundary condi- 
tions of equation (2.3.141). Now, J(t Q ) can be determined, since all com- 
ponents of its integrand have been determined. Further, the limits for the 
integral are from t Q to T, so J can be evaluated at the same time the adjoint 
equations are integrated, since t 0 varies and T is fixed. It is also conven- 
ient to tabulate 7~' It) S T (t> as the integration progresses since 

this information is needed during the execution of the guidance scheme. 

The remaining parameters in equation (2.3.152), J '(t,) U B ) , can 

also be determined while the integration of the adjoint equation is being per- 
formed. However, if the discrete-time approach is taken, these parameters 
need only be known at certain epochs of the flight. If the time t Q is taken 
as a sample point, the following mechanization can be implemented for the dis- 
crete-time case during a time interval. 
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In the previous mechanization, the state is sampled at predetermined epochs. 

The sampled state is then compared to the nominal state for the same epoch, 
and a state deviation for that epoch is determined. The state deviation then 
undergoes two multiplications. The first matrix is precalculated for the par- 
ticular sample epoch, but the second matrix is a sequence of stored matrices 
that are fed to the guidance system at such a rapid rate that it makes the 
matrix look as though it is a time-varying matrix. The remainder of the scheme 
is straightforward and will not be discussed further here. 

The mechanization of the continuous-time approach is shown below. 



The main distinctions of this type of mechanization are: (1) the precomputed 

multiplier matrix is stored as a sequence of matrices for the entire flight, 
whereas in the discrete-time case, part of this multiplier was known only at 
sampling epochs; and (2) the state variables must be known as a continuous 
function of time (or as close as possible) for the continuous-time case; the 
state was only needed at sampling epochs in the discrete-time case. 

2. 3. 2. 2. 2 Optimum Linear Control . Optimum linear control, as discussed 
here, is a method of determining control deviations by a simple multiplication 
of a precomputed matrix by the state deviation. The technique is analogous to 
a discrete-time sampled-data system in that it divides the control interval 
into a number of smaller increments and determines the control for each incre- 
ment as a separate problem. The control that is selected must then extremize 
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the performance index -which represents the cumulated performance over the 
entire trajectory. 

Since the previous control scheme and the one suggested in the preceding 
paragraph are simply different implementations of the same idea, it is expected 
that the schemes themselves are similar. First, it is noted that both schemes 
determine the control deviation by the multiplication of the state deviation 
and a precomputed matrix. Second, as it will be seen, both schemes also use 
the same type of performance index, i.e., a positive definite quadratic form. 
Thus, the primary differences exists due to variations in the methods used to 
compute the multiplier matrix in the two schemes; the results should agree, 
at least to the first-order, with those provided by the previous formulation. 

It is recalled that the multiplier matrix in the previous scheme was deter- 
mined by the multiplication of the matrix of functions (-which were inte- 
grated backwards in time from some terminal condition) and other parameters 
of the system. In contrast, the multiplier matrix in this section is deter- 
mined by a dynamic programming approach. That is, the multiplier matrix is 
determined for the last step first. ’ Then it is determined for the next to 
the last step based on extremizing the performance index for the last two 
steps. The process continues in a step by step decision process to the epoch 
of flight initiation. The result, as in the previous scheme, is a series of 
multiplier matrices for each increment. With these introductory remarks in 
mind, the optimum linear control derivation will now be presented. 

Consider the following nominal trajectory which has been divided into N 
control increments 



and consider a general performance index which accounts for state and control 
deviation for the increments from n to N 


J 




= Z [ % T U) q* XU) + U T u-0 y- uU-')\ 

4-r> 


( 2 . 3 . 153 ) 


where Qj_ and d i are arbitrarily weighting matrices determined by the control 
engineer. Finally, define the minimum loss associated with n-*-l as 
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(2.3.154) 


all u Ja 


That is, %_ n +x i.s the value that Jjj-n+l takes for the sequence of controls 
from t n _x to tjj_x that minimizes Jjj-n+i • Any other sequence of controls 
during the same time interval will produce a loss function that is greater 
than or equal to Iu_n+1 Now, if the minimum loss from t n+ ^ to tjj is L, n , 
then the minimal loss from t n to tjj is N ~ 


X 


N-n+t 


- T(h) XC»>)+ U T c^-i) UC"-') 4 X^.„J(2.3. 155) 


This statement is obvious in that it states that if an optimum control policy 
from time t n+ -, to t^, i.e.,£ii (n)... H(N-j) 3, has been determined, then the 
minimum loss for the interval [t n , tjjJ is determined by the choice of u(n-l) . 
Of all the possible choices for u(n-l), there is one that will minimize equa- 
tion (2.3.153) and the value of Ijj-n+l is thus determined. 

Now consider the loss for the last step 

X, = f UC/v-d J ( 2 . 3 . 156 ) 


But the state vector and control vectors of one time can be related to the 
next by the discrete-time form of equation (2.3.119) , i.e.. 


xch) = x c / v -/; * & N - f u(n~o 


(2.3.157) 


where 0 N _-^ — exp (At s ) = the state transition matrix 

a„_, -- f\txp ACts-V] o" 

o 

Thus, substituting equation (2. 3. 15 7) into equation (2.3.156), the loss for the 
last increment is found to be 

I, -- ^ qJ[ 0 n -, X(#-0 + A/v- 1 1 <C«-i\ ( 2 . 3 . 158 ) 
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Now, expanding and collecting terms yields 

I. 


' / > 


(2,3.159) 


so that the minimum of J^> which is I-j_, can be found by determining the solu- 
tion to 


a j, 

^ U r (A/-l) 


= o 


or 

2 A T ^, <?„ X(A/-0 * 2 [ + Y N J U(N~0 = O (2.3.160) 

Solving for the optimum control for the last step yields 

U(N~>) -- -[a*-,?***-, * y j 1 X(*-/) (2.3.161) 

or 

Cf(A/~0 - aj X(M-l) (2.3.162) 


where is defined from equation (2.3.166). 

The optimum control for the last step is thus seen to be a linear combin- 
ation of the state deviation of the (N-l)™ 1 step, /if the previous state 
deviation must be estimated, it can be shown that the optimum estimate of 
x(N-l) suffices for x(N-l)._/ 

Equations (2.3.161) and (2.3.167) can now be used to determine the minimum 
loss for the last step. 

z , --{ X r (»-0 <?» <*-') + A*-/ Qas <P*- t X(#~0 

y x^as-Og^a 7 *,-, 9/v 4/vv * *] a, 7 X 

- x t 6 v-/)[ <f> r N i q N 4 > n _, y 2 «, Qh 

(2.3.3.63) 


i- a. 


* X r (AJ-0 p, X(hl-I) 
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It is thus seen that the minimum loss for the last step is a quadratic func- 
tion of the state deviation at time N-l (the beginning of the last step) . 

The optimum control for the last two steps can be found by employing 
equation (2.3.150). 

It = U(N-Z){ +£/ T (*-z) * I, J ( 2 . 3 . 164 ) 


but 


so 


or 


*2 * 


But 


J/ = yJ(N-t) 

Z z *U(N-Z) f Z r(A/ - l) <?*/-! x. (#-') + U T (#-2) 

*- (2.3.165) 

f X r (A /-0 P, f 

UC^){^~0[q^ + Py'j X(A/~/) + U r (N- 2 ) U ( A/-z) j- ( 2 . 3 . 166 ) 
XCA /-0 = ^_ z X(A/-Z)z UCA/-Z) ( 2 . 3 . 167 ) 


Therefore, as before, the optimal u(N-2) can be from the solution of 

a j 2 


d U r (AJ-Z) 


- o 


[q„-,+ Pj] $„_ z * (/V - Z) ( 2 . 3 . 168 ) 


Finally, 
c/(a/-z) - 


( Q/v-/ * /C> ') -2 * r^,] * P?) x.CM'2) 


or 


uc/z-z) - al x(n-z) 


( 2 . 3 . 169 ) 
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and the minimum loss can be found from equation ( 2 . 3 . 166 ) as 

I 2 = X T Ct/-t) J + j X iN-/) *X r (//-2) q z C( T Z (2.3.170) 

Now, recalling that 

Xi/J-0 - <P M - Z X(V-Z) * 4y- 2 

= A «-z Q l X("- 2 ) (2.3.171) 

reduces equation (2.3.170) to 

$l_ 2 * X 7 («-2)« 2 K-i] [<?#-, + P *][0#-2 Jf^-2V4^^(^2)] 

+ X 7 C N- z) <l z tr*-, #2 X ( /V- 2) 

■ { ^-2 [Q#-, i ‘ p >] 2 F , *1 } 1 Ca/- 2) (2.3.172) 

where 

F > = *l- 2 ( Q#-t + P,) */v-z -/ 2T^ X (2.3.173) 


Thus, writing equation (2.3.172) in the form of equation (2.3.163) 

T z - X \m-z) P z X_(N-Z) (2.3.174) 

where 

P Z - 0Z-Z (Q«-, "?) <S> N _ Z -Q z F, 4 z r (2.3.175) 

This procedure can be extended to the n^* 1 step as follows: 

CfCn) = aT Ai-y> (2.3.176a) 
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' ~ F *'n ^ T »-/[Qn + p H-n ] ^-/ 

^/i r * ^] 

^ N-nhl '- $»-* [ ^ * ^-/l] 4''-’ ~ Q H-'>+/ Fa/ -» 

P<? - o 


(2.3.176b) 

(2.3.176c) 

<^ 2 ‘ 3 * 176d) 

(2.3.176e) 


The use of these equations should be apparent from the previous discussion. 
The feedback matrix, a? , can be found by starting with n=N and iterating 
back to n=0. “ n 

The sequence of steps is shown below: 

1. , Compute F c 


_/ u 

P<0 ~ & A/-/ C Q /V ^y?o) A N-/ * 

2. Compute a£ 

~ ~ Po A T w~t ( Qn * ^o') Pm-/ 

— ~ a/-j Qn Av-/ * <^v ^/V-/ Qh &N-/ 

3. Control for time t^_^_ to tjj is 

u(n-i) r a T t x(/v-/) 

4. Compute 

p, -- &.,[<?+/£] - <*, ^0 < 


5. Compute F-j_ 


f-j - 


- A m-2 An-2 + 
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6. Compute 


< = 




7. Control is 


UiAt-Z) - Cc± XC//- 2) 


The general flow chart for this process is 


N -/? 


N-n+/ 


A N~/ (Qn * pN-n) <f. 




The actual mechanisation of this optimal guidance control scheme would consist 
of the storage of all of the elements for the "a" matrices from time t]_ to t„. 
The control deviation can then be established by a multiplication of the state 
deviation for the particular time of interest by the current feedback matrix. 


An example of the application of optimum linear control follows: Let 

£ r and5r be the position. and velocity perturbations from 3ome nominal trajec- 
tory. The deviations at time t, can be related to those at time t Q as follows 








(2.3.177) 


o> 1 
1-^ 

1 


1 

0/ ' 02 


ifl 

CO 

1 

■h 

0z 



1 


8 r G 
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where £ Vq is some velocity correction that can be made at time t Q . This 
equation can be thought of as the discrete-time form of the general solution 
to the linear system 


X(/) - + A a U(0) 


(2.3.178) 


A loss function for this process can be defined as the square of the position 
deviation at time t]_, i.e.. 




r 1 n 

I \ o 


sn 

1 

r • *| 




7 , -- &n • = 

6 Jj J 

0 \ o 


s'n 


(2.3.179) 


Under this particular choice of a loss function, no penalty is placed on the 
fuel needed in order to perform the control correction. In "real life" prob- 
lems, this choice of a loss function would be less than ideal since the loss 
function in effect forces the position deviation to be minimum at any fuel 
cost. Ideally, it is desirable to use a loss function that not only weighs 
the position deviation, but also the control deviation, since the nominal tra- 
jectory is optimal. 

Using the standard notation, the following definitions can be made from 
equations (2,3.177) and (2.3.178). 



The loss function, equation (2. 3. 179 ) > may also be written as 


Jz -- ZjQ'h 


(2.3.180) 


118 


so 


Q - 


I 

o 


o 

o 


r, = o p D = <0 


Thus, the previous results of this section, in conjunction with consideration 
of the problem of one interval, define the control as follows: 


££(0) = <*, r X(o) 




Ao 9,0o 


(2.3.181a) 

(2.3.181b) 


F B - 


-- q,a 0 .#] 


O * D 

u « 


a 


= 0\0 Z (2.3.181c) 




*: 




i , o 

I — i — i 
O • o 

L i _l L i 


\<t>, \ 0 M ' 
0~ j" i/ 


~\^I <**] [0 e r jA] 


0, I 02 

. _l 

o ' o 


(2.3.182b) 


The optimum linear control is therefore 


«« ■«'(!!) 


or 


6 Vo = - <0/ [<A ! 02. ] 


5 If 

8 r e 


(2.3.183) 


(2.3.184) 


s Vo -- - 0/0, *{o - si 


(2.3.185) 
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A simple verification 
position deviation, <£ r-^. 

s n - <P, & £ 

= & 5£ 0 
= O 

It should be noted that the results to this over-simplified problem are iden- 
tical to the "Fixed Time of Arrival" scheme that is presented in Reference 
3.25. The reason for the similarity is that only one control increment was 
analyzed and that the only performance criteria for the mission was the posi- 
tion deviation at the end of the increment. The answer to the previous prob- 
lem is not so obvious when a more complicated loss function is considered or 
when a "loss" is associated with the control deviation. 


of this answer can be made by determining the 
Substitution into equation (2.3.175) yields 


* 8£o i- < 51 /* 

f ^2 [“ 0* ^ ~ * ts j (2.3.186) 
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3.0 RECOMMENDED PROCEDURES 

The material presented in the previous discussions is applicable to a 
series of boost vehicle guidance problems. However, since the guidance 
equations are but a small part of the total guidance loop and since the 
selection of the guidance equations should be based upon consideration of 
the required performance from the total loop, the reliability and cost of 
the system, the mission success criteria, etc., any attempt at this point 
to select a specific set of guidance equations for general application will 
fail. This observation is strengthened when it is realized that the primary 
differences in the results of these applications as observed in terms of the 
mission is in the flexibility derived. Thus, since no set of guidance 
equations appears to have a clearly defined advantage for all applications 
prior to consideration of the total guidance loop, a definitive decision 
will be deferred until the monograph on guidance system synthesis is pre- 
pared. However, general preferences for the guidance of boost vehicles 
of the present and near future can be presented. 

The simplest of the schemes (perturbation guidance) is applicable in 
varying degrees of sophistication to the boost of ballistic missiles, small 
scientific payloads, and manned satellites (where the trajectories are well 
defined before launch and where retargeting has been ruled out as a possible 
event). In contrast, while the more complex mechanizations of the explicit 
and adaptive guidance approaches do not preclude application to the same 
missions, the cost of the system, resulting from more elaborate instrumenta- 
tion said a larger guidance computer, will generally preclude application to 
those missions which do not require flexibility. Thus, the adaptive mode 
of guidance will probably be restricted, for the time being, to the more 
demanding missions such as manned lunar and planetary escape trajectories, 
and extraterrestrial soft landing. This application is justified for these 
mission phases due to the complex manner that the mission phases are tied 
together and due to the variety of abort and mission redefinition possi- 
bilities which exist (to derive usefulness from a partial success). 

Finally, in the case of the adaptive type of guidance, the two forms 
would appear to have the same type of computational requirements and appear 
to impose the same order of mechanizations!, complexity. However, since the 
iterative path adaptive guidance mode of section 2.1 requires less targeting 
before the vehicle is launched, this approach is preferred for the more 
sophisticated applications over the guidance polynomial approach of section 
2 . 2 . 


Mechanizations of the material presented in both sections 2.1 and 2.3 
are illustrated in the respective sections. Thus, reference is made to these 
sections for such information. These discussions show the unique capabilities 
of the guidance equations and illustrate the degree of optimization attained 
in the process. 
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APPENDIX A 


Formulation of the Variable Mass Vehicle Problem 


This appendix will present a straightforward formulation of the variable 
mass vehicle and then show why a different selection of the state variables 
is desirable in order to avoid cross coupling the state variables in the 
solution. A conventional formulation would assign the state variables as 


follows: 


= 


K.- 


B xhaust V&jLoary 


z z = x s 


u; 




+ «? =/ 


^ * 


where 



~>n(t) 

7n it) 



'*h Ct) 
•Jnit) 
7n(t) 


U* 


U 4 


(i) - vehicle mass 


(A.l) 


Z, j 5 Xj - position of vehicle with respect to target 
* z s - » x 6 = velocity of vehicle with respect to target 


In addition, the following state and control definitions are made 

•»<t) = x 7 
**U) = 

So the state equations become 


A ~ *v 

= *r 

*3 ~ ^ 



(A.2) 


(A.3) 
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X S ~ Hs 

* . 

* ' ^5 

6 U 7 

where the admissible control is 

° ~ ^(t) ^ D 

This formulation seems reasonable since the controls a., , u.^ , u_ 3 determine 
the direction of thrust and the control u ¥ determines the magnitude. How- 
ever, it will be shown to have an over complicated solution once the Maximum 
Principle is applied. This fact is seen in the following analysis. 

It is required to determine the control vector components a. , u. x , u 3 
and Uy such that the relative position and velocity between the vehicle and 
some target are reduced to zero while x 7 Ct = T ) is a maximum. Applying 
Pontry agin S Maximum Principle, it is desired to maximize 5 r = c - X 7 where 
C 7 = / and the remaining components of c are to be determined. For the 
system of equations, (A. 3), the function Hi % } P , u. , -L ) becomes 


<4 


Mix.EyU.t) = p,X^ p z X s tp 3 x b r 

A i\ 

The system T , thus becomes 


U >P« +e *Z.Ps +“jP 6 ~-7T p > 
ye 


(AA) 


A - P 3 ~ 0 


£,= -p> 
Ps~ 'Pz 

A = 


4= ^ ^[^v***^*^] 


(A.5) 


Now, according to Pontryagin’s Maximum Principle the optimum control vector 
is determined by minimizing H ( % , B , «• , i ) . However, equation (A. 4) 
and (A. 5) are somewhat "over” complicated since the determination of the 
optimum control vector components is dependent on X 7 and P 7 , 
n.nfl p , in turn, is dependent on the control vector components and 

% 7 7 . This "over" complication of equations (A . k and A.5) resulted 

from the definition of "X 7 as the mass of the vehicle. It is seen that 
the process of state variables for the system which involve products of 
system states and control vector components should be avoided wherever possible. 
Fortunately, it is possible to do so in the present problem if an alternate 
control vector component u. is defined such that 
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where the admissible control is 

Oi U r (t)± D 


This formulation seems reasonable since the controls u, , u.^ , u_ 3 determine 
the direction of thrust and the control u ¥ determines the magnitude. How- 
ever, it will be shown to have an over complicated solution once the Maximum 
Principle is applied. This fact is seen in the following analysis. 


It is required to determine the control vector components u, , u. x > u 3 
and such that the relative position and velocity between the vehicle and 

some target are reduced to zero while % 7 (<= T ) is a maximum. Applying 
Pontryagin Maximum Principle, it is desired to maximize 5 r = c where 

C 7 z / and the remaining components of c are to be determined. For the 
system of equations, (A. 3), the function H ( £ ,P , u ,t ) becomes 


"fz, = P,x« + PzXs+PjXt,* Ve V 

The system P , thus becomes 


u 'P*, +i <iPs '■“Sr-rrP' 

re 


(A. 4) 


P, m Pz “ P 3 ~ ° 


4 = ~p 

Ps “ 'P* 
a = -Ps 

4 = * 


<V v r 

"~t [4 + *** pf + 


(a.5) 


Now, according to Pontryagin' s Maximum Principle the optimum control vector 
is determined by minimizing P i ?< F j , i) . However, equation (A. 4) 
and (A. 5 ) are somewhat "over" complicated since the determination of the 
optimum control vector components Is dependent on X. ? and P 7 , 
and p , In turn, is dependent on the control vector components and 
"X . 7 7 . This "over" complication of equations (A. 4 and A.5) resulted 

from the definition of % 7 as the mass of the vehicle. It Is seen that 
the process of state variables for the system which involve products of 
system states and control vector components should be avoided wherever possible. 
Fortunately, it is possible to do so in the present problem if an alternate 
control vector component u is defined such that 
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*1- x s 
*3= *6 


*5 ' U S U Z 
*4 = V e “s U s 
x 7 = K* u s 


(A.6) 


The control vector component u ff is introduced such that in equation (A. 3) 

is essentially absorbed. It is important to note that ih equation 

(A. 6) is no longer the mass of the vehicle since u s is no longer the mass 
flow rate. Actually 7 u s is defined as 


«s U) = " 


^ U ) 
7* It) 


(A.7) 


Thus, 


z 7 (r) = x 7 (o) 


dt 



d~m 

-m 


Z 7 M = z 7 (0) - V e Jly a ™U)\ o 


"2>? (Of 
yr?('r) 


(A.8) 


where X(o) is arbitrarily selected as zero without loss of generality. 
Thus, for minimum mass expenditure it is desired to minimize t) or, 

equivalently, it is desired to minimize S T * l - x 7 (t) . 

In terms of U. g , the function H( x ,p , tt J) and the P system become 


1 ) = P> *v * P, + Pj x 6 


" s \u, p« + v 2 p r +u 3 p 6 +p 7 ] 


(A.9) 


P, - p z = p j = o 

Ph =-p, 

Ps = (A.IO) 
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Now, since " O * py (£) ~ p y (?) - ~~f 


> 


therefore. 


Mi*, e, <*,*)* r.Zf+Pt x s*PjX, * v c u s *u, />,<■ u 3 ? t -/] (A<11) 

A comparison of these equations with equations (A. 4) and (A. 5 ) shows the ad- 
vantage gained in the use of ol s as a control vector component and the 
alternate definition of X 7 

The previous formulation was part of a more general study on the applica- 
tion of Pontryagin * s Maximum Principle to optimum control of a variable mass 
space vehicle performed at HAA by D. R. Grier. 
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